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COCHISE  Studies  I: 
Fluid  Dynamical  and 
Infrared  Spectral  Analyses 
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1.  INTRODUCTION 

The  COCHISE  facility  consists  of  a cryogenically  cooled  test  chamber,  provid- 
ing operating  conditions  of  P ~ 10”^  torr  and  T = 40  - 100°K,  that  was  designed 
for  the  study  of  chemiluminescent  reaction  kinetics  under  essentially  "single"  col- 
lision conditions.  Reactant  gases  are  introduced  into  the  test  chamber  through 
two  sets  of  diametrically  opposed  inlet  tubes.  The  gas  entering  through  one  of 
these  sets  can  be  passed  through  a ^ -wave  discharge  to  produce  free  radicals  or 
electronically  excited  species.  The  concentrations  of  such  species  can  be  moni- 
tored via  absorption  spectroscopy.  The  primary  diagnostic  of  the  system  is  an 
infrared  spectrometer  used  to  monitor  the  spectrally  resolved  radiation  arising 
from  the  reaction  under  study.  The  resulting  radiation  histories  can  then  be  inter- 
preted in  terms  of  reaction  quantum  efficiencies  and  the  excited  state  distributions 
of  reaction  products. 

The  research  conducted  during  the  reporting  period  is  broken  into  two  parts: 

The  first,  described  in  Section  2,  involved  the  development  of  a detailed  fluid 

dynamical  analysis  of  the  flow  field  and  operating  conditions  of  the  test  chamber. 

The  second  part,  described  in  Section  3,  includes  a brief  examination  of  a proposed 

2 

absorption  technique  for  measuring  the  concentration  of  N(  D)  atoms,  and  a review 
and  interpretation  of  preliminary  measurements  performed  on  the  system. 

(Received  for  publication  6 December  1977) 
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2.  COCHISE  FLUID  DYNAMICS 


2.1  System  Definition 

A schematic  diagram  of  the  COCHISE  facility  is  shown  in  Figure  1.  Two  rows 
of  diametrically  opposed  inlet  jets  carry  the  relevant  chemical  species  into  the 
tank.  The  chemical  kinetics  of  the  species  are  monitored  in  the  region  within 
2 cm  of  the  tsmk  centerline.  The  gas  then  flows  toward  the  cryogenically  cooled 
walls  which  provide  the  necessary  pumping  capacity.  Nominally,  wall  tempera- 
tures of  20°K  to  40°K  are  expected  to  provide  chamber  pressures  in  the  range  of 
10  ^ torr. 

The  inlet  gas  flow  is  controlled  with  a needle  valve.  Approximately  25  cm 
upstream  of  the  tank  inlet,  the  gas  pressure  (p.),  temperature  (T.),  and  mass  flow 
per  inlet  tube  (M)  are  measured.  With  this  information,  along  with  the  wall  tem- 
perature and  the  tank  geometry,  the  temperature,  pressure,  and  density  of  the 
gas  on  the  tank  centerline  can  be  deduced.  In  addition,  a general  flowfield  visual- 
ization can  be  developed. 

A schematic  representation  of  the  flowfield  is  shown  in  Figure  2.  Top  and  end 
views  and  a basic  recirculation  cell  for  each  jet  are  illustrated.  Each  cell  is 
divided  into  four  regions.  Region  1 is  the  inlet  pipe  and  region  2 is  the  expansion 
of  the  inlet  flow  into  a relative  vacuum.  Region  3 is  an  axisymmetric  stagnation 
point.  The  flow  in  region  2 is  supersonic  and  must  become  subsonic  to  match  the 
stagnation  point  flow.  The  structure  surrounding  region  2 is  meant  to  depict  a 

barrel  shock  and  Mach  disc,  such  as  that  surrounding  a high  altitude  rocket 
1-3  4 5 

plume  or  free  jet  expansion.  ’ These  shocks  transform  a high  velocity,  low 
pressure  flow  into  a low  velocity  flow  of  higher  pressure.  Region  4 is  also  sub- 
sonic, but  the  appropriate  model  is  not  yet  apparent. 

The  total  temperature  at  the  stagnation  point  in  region  3 is  the  total  tempera- 
ture entering  the  tank  through  the  inlet  jets.  The  stagnation  pressure  (p^)  in  reg- 
ion 3 and  the  location  of  the  shocks  are  interrelated.  We  must  determine  a value 
of  Pq^  that  is  compatible  with  the  cryogenic  pumping  at  the  walls,  then  determine 
the  location  of  the  shocks  from  p^^.  Hence,  p^^  is  determined  only  from  the  wall 
capture  coefficient  C.  These  coefficients  are  determined  empirically^  under  the 

1.  Albini,  F.A.  (1965)  AIAA  Journal  3;1535. 

2.  Hubbard,  E.  W.  (1966)  AIAA  Journal  4;1877. 

3.  Boynton,  F.  P.  (1967)  AIAA  Journal  5:1703. 

4.  Ashkenas,  H. , and  Sherman,  F.  S.  (1966)  The  structure  and  utilization  of 

supersonic  free  jets  in  low  density  wind  tunnels,  in  Rarefied  Gas  E^amics. 
Fourth  Symposium,  Vol.  11.  Academic  Press,  New  York,  pp.  84-105. 

5.  Bier,  K. , and  Schmidt,  B.  (1961)  Zeitschrift  Agnew  Physik  13:495. 

6.  Dawson,  J.  P.  (1966)  J.  Spacecraft  3!218. 
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Figure  2.  Schematic  Diagram  of  COCHISE  Flowfield 

assumption  that  the  gas  in  the  tank  is  static  with  temperature  Tq^,  pressure  Pq^, 
and  density  To  develop  a more  detailed  model  for  region  4 and  apply  capture 

coefficients  from  Ref.  6 would  be  Inconsistent  and  would  probably  introduce  extran 
i eous  errors.  We  shall  adopt  the  model  for  region  4 that  is  consistent  with  the 

manner  in  which  the  capture  coefficients  were  measured;  that  is,  the  gas  is  taken 
to  be  static  with  a uniform  pressure,  temperature,  and  density. 

A more  detailed  model  for  region  4 is  developed  in  iVppendix  A where  the  gas 
is  shown  to  be  a viscous  continuum  flow  described  by  the  Navier-Stokes  equations 
with  "slip  flow"  at  the  wall.  However,  this  model  was  abandoned  when  it  was 
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realized  that  the  capture  coefficients  were  not  known  accurately  enough  to  justify 
such  a detailed  analysis  of  the  flowfield. 


I 


2.2  Inlet  Tubes 

The  gas  temperature,  pressure,  and  mass  flow  in  the  inlet  tubes  are  meas- 
ured approximately  25  cm  before  entering  the  tank.  Let  us  denote  these  conditions 
as  Ti,  Pj,  and  M,  respectively.  The  gas  density  corresponding  to  p.  and  T.  is 
denoted  by  and  the  gas  velocity  is  u^, 

u.  =4M/p.trd^  , (1) 

where  d is  the  tube  diameter  (1.  1 cm). 

We  wish  to  determine  the  conditions  at  the  end  of  the  inlet  tube,  where  the  gas 

enters  the  tank.  Let  us  denote  these  conditions  with  a subscript  unity.  The  pres- 

7 

sure  drop  corresponding  to  a laminar,  incompressible  pipe  flow  is 


XL  2 

Pi-Pl=■2dPi^  ' 

where  L is  the  length  of  the  pipe  (25  cm),  X is  a coefficient  of  friction, 
X = 64/Re^  , 

and  Re^  is  a Reynolds  number, 

Rej=pju.d/pj  . 

From  the  above  expressions,  it  follows  that 
p,  128Mp.afL 


(2) 


(3) 


(4) 


(5) 


where  a.  is  the  inlet  speed  of  sound. 

The  flow  Reynolds  number,  Re^,  may  be  expressed  as 

Rej  = 4M/»dMj  (6) 


7,  Schlichting,  H.  (1955)  Boundary  Layer  Theory,  McGraw  Hill. 
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and  the  Mach  number,  Mj,  may  be  written 

M.  = 41VIa./7irp.d^  . (7) 

Equations  (5)  through  (7)  impose  three  constraints  on  the  combinations  of  M 
and  Pj  that  are  compatible  with  the  inlet  pipe.  We  require  Pj^  > 0,  < 1,  and 

Re^  < 2300.  The  Reynolds  number  constraint  is  imposed  to  prevent  turbulent  pipe 
flow.  These  constraints  are  illustrated  in  Figure  3 for  Tj  = 100®K.  Clearly  we 
should  keep  our  pressure  losses  to  less  than  10  percent  to  avoid  significant  viscous 
heating  of  the  inlet  gas.  Maintaining  Pj^  > 0.  9 p.,  M^  < 0.  25,  and  Re^  < 2300,  we 
can  be  sure  that  Pj^  and  Tj^  are  within  10  percent  of  our  measurements  of  p^  and 
Tj.  Figure  4 illustrates  the  conditions  imder  which  these  constraints  are  satisfied. 
The  calculations  were  carried  out  for  Tj^  = 300°K,  100°K,  and  60°K,  corresponding 

to  a.  = 30,  000  cm/sec,  17,  000  cm/sec,  and  13,  000  cm/sec,  respectively,  and  for 
^ -4  -4  -5 

= 2 X 10  ,10  , and  6 X 10  gm/cm-sec,  respectively.  If  a pressure  of  5 to 

10  torr  can  be  maintained  in  the  inlet  pipe,  the  conditions  measured  in  the  inlet 
tube  are  equivalent  to  those  under  which  the  gas  enters  the  tank. 

In  the  event  that  the  5 to  10  torr  pressure  cannot  be  maintained,  the  viscous 
losses  will  heat  the  gas.  This  heating  can  be  offset  by  maintaining  the  walls  of  the 
inlet  pipes  at  temperature  Tj.  The  question  is,  how  hot  will  the  gas  become  before 
the  dissipation  will  be  balanced  by  heat  conduction  to  the  wall?  Let  us  seek  the 
final  state  where  conduction  does  balance  dissipation.  We  write 


(8) 


or 


(9) 


To  determine  the  constant  in  the  above  relationship,  we  acknowledge  that  viscosity 
cannot  heat  the  gas  above  the  total  temperature.  Therefore, 


AT 

T. 


(10) 


Since  Mj  < 1,  the  maximum  possible  temperature  rise  in  the  inlet  tube  is  20  per- 
cent for  7 = 1.4. 
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2.3  Free  Expansion  Region 


The  flow  within  the  barrel  shock  and  Mach  disc  expands  as  though  there  were 
no  back  pressure.  The  gas  reaches  a limiting  velocity  in  which  all  of  the  thermal 
energy  of  the  flow  is  converted  into  kinetic  energy, 


(11) 


Since  the  inlet  Mach  number  is  subsonic,  it  follows  that  « 2/(7  - 1)  and  the 
stagnation  temperature  in  the  free  jet  is  approximately  T..  Hence, 


= T. 


(12) 


from  which  we  obtain 


'^2  = 

where  a.  is  the  speed  of  sound  corresponding  to  T., 

*•  '34 

The  gas  density  in  the  free  expansion  is  approximately  ’ given  by 


P2 


Af(e)M 


(14) 


where  A is  a constant,  r2  is  the  spherical  radius  from  the  exit  plane  of  the  inlet 
tube,  0 is  measured  from  the  jet  centerline,  and  [(g)  is  obtained  from  numerical 

3 

calculations  of  high  altitude  rocket  plumes.  A good  approximation  for  [(g)  is 


He)  = 


(^) 

X ^ 


2/(7-1) 


(15) 


where  g is  the  maximum  turning  angle  for  the  expansion  of  an  isentropic  gas  to 
a vacuum  (130.  6°  for  7=1.  4). 

Conservation  of  mass  requires  that  A be  determined  by 


-1 


[(g)  sin  g do 


(16) 
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Since  most  of  the  mass  flux  occurs  near  0 = 0,  we  approximate  the  sine  of  0 by 

(17) 

and  determine  A to  be 

A = <'y  + , . (18) 

8(7  - 1)  0“, 

00 

The  expressions  for  the  gas  velocity  Ug  and  the  gas  density  pg  are  valid  within 
the  shock  structure.  The  location  of  the  shock  structure  must  be  determined  as  a 
fimction  of  the  stagnation  pressure  on  the  tank  centerline. 

2.4  Stagnation  Point  - Region  3 

The  diametrically  opposed  inlet  jets  will  interact  in  a plane  of  symmetry  pass- 
ing through  the  tank  centerline.  This  was  illustrated  in  Figure  2.  Symmetry  about 
the  jet  centerline  indicates  that  the  flow  may  be  modeled  as  an  inviscid,  axisym- 
metric,  stagnation  point.  Furthermore,  in  passing  from  region  2 to  region  3,  the 
flow  on  the  jet  centerline  must  traverse  a plane  normal  shock  wave.  Hence,  the 
flow  in  the  stagnation  region  is  subsonic  and  may  be  treated  as  incompressible. 

In  addition  to  the  incompressible  assumption,  we  assume  that  the  stagnation  point 
solution  is  equivalent  to  that  for  which  the  flow  is  initially  normal  to  the  plane  of 
symmetry.  This  approximation  is  strictly  valid  only  on  the  jet  centerline,  where 
regions  2 and  3 will  be  matched. 


The  incompressible,  potential  flow  solution  for  the  axisymmetric  stagnation 

7 

point  as  described  above  is  given  by 


Figure  5.  Stagnation  Flow  for  Region  3 


Denoting  the  location  of  the  shock  by  z = Z^,  we  note  that 


m = ■ 


2Z. 


(22) 


where  Ug  (Z^)  is  the  post  shock  value  of  Ug.  For  strong  shocks  (which  is  compat- 
ible with  the  assumption  that  Ug  is  a limiting  velocity),  we  write 


(23) 


Therefore, 


m = (y  - 1)  U2/2  (y  + 1)  Z^ 


(24) 


where  Z^  is  unlinown  and  must  be  determined  through  knowledge  of  the  cryogenic 


pumping  at  the  chamber  walls. 
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Denoting  the  distance  from  the  exit  plane  of  the  inlet  tube  to  the  Mach  disc  by 
Rg(0),  we  note  that  Rg(0  = 0)  is  related  to  2^  by 

Rg(0)  = X - Zg  . (25) 

where  X was  defined  in  Figure  1.  The  total  pressure  in  the  stagnation  region 
becomes 


= 


(y  + 1)  JTMUg 
8(y  - 1)  (X  - Zy 

00  9 


(26) 


These  results  for  the  total  pressure  are  illustrated  in  Figure  6.  For  M = 0.  03  gm/ 

_3 

sec  per  inlet  tube,  pressures  of  the  order  of  10  torr  can  be  obtained.  It  must 
be  remembered  that  the  location  of  the  shock  is  a function  of  the  pressure,  p-  , 

°3 

and  that  pressure  is  a function  of  the  cryogenic  pumping.  However,  Z is  an 


Figure  6.  Stagnation  Pressure  on  Tank  Centerline 
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important  parameter  in  the  analysis,  because  we  must  keep  the  shock  a couple  of 
centimeters  away  from  the  tank  centerline  to  ensure  that  the  optical  monitor  of  the 
experiment  does  not  see  the  underexpanded  flow  of  region  2.  Physically,  can- 
not become  zero  since  continuity  requires  that  some  minimum  value  of  must 
exist.  Raising  the  background  pressure  could  increase  above  its  minimum 
value,  but  lowering  the  background  pressure  would  result  in  Z^  monotonically 
approaching  a minimum  value  to  be  determined  by  a Newtonian  pressure  balance. 
The  Newtonian  balance  is  illustrated  in  Figyre  7.  The  inlet  tube  acts  as  a spherical 
source  that  is  reflected  from  the  plane  of  symmetry.  The  component  of  velocity 
normal  to  the  plane  of  symmetry  is  lost  to  pressure,  but  U2  sin  $ is  conserved 
across  the  shock  and  moves  the  high  pressure,  low  velocity  gas  lateral  to  the  jet 
centerline.  Mass  balance  suggests 


*’2^2’^ ''2  ' ^®P2^^'^2  21^26 


6 = X/12  = 1.  7 cm 


(27) 


(28) 


Thus,  the  minimum  value  of  Z is  1.  7 cm  and 

s 

Zg  > 1.7  cm  , (29) 

To  determine  the  static  pressure  along  the  tank  centerline,  we  require  knowl- 
edge of  the  temperature  and  density.  The  total  temperature  is  conserved  across 
the  shock  wave. 


The  total  pressure  and  temperature  yield  the  total  density. 


y(y  + 1)  IT  M 

Pn ^ 15 3“ 

°3  4(7  - 1)^  u,©*^  (X  - ZY 

i 00  S 

and  the  expression  for  the  static  pressure  becomes 

Pg  y(y  - l)(rg  + 4 z^) 

P03  ■ ' ■ “4(7  + 1)2  Zf 


(30) 


(31) 


(32) 


where  we  have  assumed  that  pg  is  approximately  p^  . The  results  for  the  static 
pressure  are  illustrated  in  Figure  8.  The  closer  the  shock  wave  is  to  the  tank 
centerline,  the  faster  the  pressure  decays  along  the  centerline.  The  fastest  decay 
occurs  for  the  Newtonian  value  of  Z^.  For  = 1.  7 cm,  the  pressure  would  decay 
to  a minimum  of  60  percent  of  the  stagnation  value  at  the  edge  of  the  recirculation 
zone.  However,  this  position  on  the  tank  centerline  is  also  a stagnation  point. 

The  static  pressure  will  return  to  p^  , as  indicated  by  the  dashed  lines.  Hence, 
there  is  only  a 10  percent  variation  in  pressure  along  the  tank  centerline.  We  will 
now  proceed  to  determine  p^^  from  the  cryogenic  pumping  and  then  determine  Z^ 
from  Po  . 


2.5  Cryogenic  Pumping  • Region  4 

The  relative  vacuum  within  the  chamber  is  maintained  by  cooling  the  walls  of 
the  chamber  to  the  point  where  molecules  striking  the  walls  have  a high  probability 
of  sticking  in  a condensed  phase.  The  capture  coefficient  C for  "pre-frosted"  walls 
has  been  measured  by  Dawson.  ® He  assumes  that  the  gas  within  the  chamber  is  a 
stationary  Maxwellian  gas  of  temperature  T.,  pressure  p„  , and  density  . The 

1 O3  *'  ^ O3 

maximum  mass  flux  toward  the  wall  per  unit  area  is  determined  from  kinetic  theory: 


« 
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r j - Distance  From  Jet  ({^  on  Tank  ([  (ctnl 


Figure  8,  Pressure  on  Tank  Centerline 
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/a.  . 


(33) 


The  actual  mass  flux  toward  the  wall  is  used  to  define  the  capture  coefficient  C, 


C = 


actual 


m 


(34) 


A theoretical  basis  for  determining  C is  called  the  critical  velocity  model.  It 
states  that  there  is  some  maximum  normal  velocity  with  which  a molecule  can 
strike  a pre-frosted  wall  and  stick.  In  a Maxwellisui  distribution  function,  there 
always  exist  some  fraction  of  the  molecules  that  are  capable  of  sticking.  That 
fraction  is  denoted  by  C: 


C = 1 


exp  (-Tc/Ti) 


(35) 
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The  quantity  r is  a critical  temperature  corresponding  to  the  critical  velocity, 
fi  c 

Dawson  has  measured  for  various  gases,  and  his  results  are  shown  in  Table  1. 


Table  1.  Critical  Temperature  for 
Sticking 


In  reality,  the  behaviour  of  the  gas  in  the  tank  is  not  as  idealistic  as  Dawson 
assumed.  Once  molecules  are  captured  at  the  wall,  there  is  a mean  velocity 
toward  the  wall.  Hence,  the  Maxwellian  distribution  fimction  near  the  wall  should 
include  a mean  velocity  u^,  and  a gas  temperature  T^  that  is  different  from  T.. 
This  distribution  function  would  yield  a smaller  fraction  of  molecules  with  normal 
velocity  less  than  "critical. " In  addition,  the  basic  expression  for  rh  would 
differ.  Hence,  the  capture  coefficient  would  be  different  from  those  given  by 
Dawson. 

While  it  is  "desirable"  to  give  a more  detailed  description  of  the  COCHISE 
facility,  one  does  not  know  the  capture  coefficients  accurately  enough  to  justify 
such  an  analysis.  The  most  realistic  approach  is  to  use  the  model  that  was  as- 
sumed in  the  measurements  of  C,  that  is,  that  of  a static  Maxwellian  gas.  The 
more  detailed  model  is  described  in  Appendix  A and  most  likely  will  not  be  used 
unless  there  are  some  gross  Inadequacies  in  the  static  model. 

Using  the  same  assumption  as  Dawson,  we  conserve  mass  flux.  That  is,  the 
mass  flowing  into  the  tank  through  two  opposing  jets  must  be  captured  over  the 
wall  area  2 xR  Y,  where  R is  radius  of  the  tank; 
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Equation  (37)  yields  the  stagnation  pressure  as  a function  of  mass  flow,  inlet  tem- 
perature, and  tank  geometry.  The  primary  constraint  on  p^^  is  that  it  must  be  in 
excess  of  the  gas  vapor  pressure  at  the  wall  temperature  t^. 


p > p (t  ) , 

*^Og  *^v  w ’ 


(38) 


otherwise  the  condensed  phase  would  vaporize  and  raise  the  pressure  within  the 
tank.  For  most  gases 


inp  - A - B/t 
t-v  ' w 


(39) 


where,  for  Ng,  A = 16.  9,  B = 800  K,  and  p^  is  in  torr. 

Equations  (37)  and  (39)  yield  the  initial  predictions  for  the  tank  centerline 
pressure.  Figure  9 illustrates  Pq^  as  a function  of  t^,  T^,  and  M.  The  wall 
temperature  is  assumed  to  have  no  influence  on  the  stagnation  pressure,  provided 
Pq^  > 10  p^.  For  larger  wall  temperatures,  the  pressure  is  assumed  to  increase 
and  follow  the  vapor  pressure  curve.  The  "cold  wall"  results  are  illustrated  as  a 
function  of  mass  flow  in  Figure  10.  The  maximum  wall  temperature  for  which  the 
results  are  valid  is  indicated.  There  is  also  an  upper  limit  on  M for  each  T^  to 
prevent  the  inlet  pipe  from  becoming  turbulent.  This  is  so  indicated  in  Figure  10, 
as  obtained  from  Figure  4,  Each  value  of  T.  possesses,  within  the  framework  of 
this  model,  a unique  value  of  shock  stand-off  distance,  Z^.  This  is  obtained  from 
Eqs.  (26)  and  (37).  For  T.  of  100°K  and  60°K,  this  value  of  falls  below  the 
Newtonian  limit  of  1.  7 cm.  Hence,  the  corresponding  value  of  Z^  is  corrected  to 
1.7  cm,  as  indicated  in  Figure  10.  The  additional  constraint  expressed  in  Fig- 
ure 10  is  a gross  approximation  to  the  effect  of  diffusion  from  the  wall  to  the  tank 
centerline.  If,  at  t = o,  we  "turn  on"  the  diffusion  at  the  wall,  it  will  diffuse  with 
velocity  /t  where  v is  the  kinematic  viscosity.  The  diffusion  front  will  propa- 
gate toward  the  tank  centerline  until  the  diffusion  velocity  is  just  balanced  by  the 
mean  velocity  u.  Hence, 
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A/R=Ai7rY/M  , (42) 


Diffusion  from  the  wall  will  affect  the  centerline  when  A/R  is  order  unity.  This  is 
so  indicated  in  Figure  10  and  represents  a minimum  mass  flux  under  which  the 
tank  can  operate  without  contamination  of  the  teuik  centerline  by  the  chamber  walls. 

Figure  10  represents  the  major  result  of  this  study:  tank  centerline  pressure 
as  a function  of  the  measured  input  conditions.  Upper  and  lower  bounds  on  the 
mass  flow  for  which  the  analysis  is  valid  are  also  indicated. 


2.6  Flow  Visualization 


The  preceding  analysis  has  predicted  the  tank  centerline  pressures  as  a func- 
tion of  the  mass  flow,  gas  temperature,  and  wall  temperature.  The  analysis  also 
predicts  the  location  of  the  shock  waves  with  respect  to  the  tank  centerline.  What 
remains  is  a scale  drawing  of  the  shock  structure  within  the  tank  and  a visualiza- 
tion of  how  the  gas  flows  through  the  shocks  and  to  the  walls. 

The  location  of  the  shock  from  the  exit  plane  of  the  inlet  tube  on  the  jet  center- 
line  was  determined  in  Section  2.4, 


Rs(0)  = X - 


(25) 


2 . 


To  determine  Rg(0),  we  find  the  locus  of  points  such  that  P2(0)n2  is  identical  to 


p . Thus, 


Rg(9)  = Rg(0)  yi(^ 


(43) 


Rg(5)  is  shown  in  Figures  11  and  12.  Figure  11  shows  that  the  free  jet  expansions 


actually  impinge  and  minimize  the  flow  along  the  axis  of  the  tank,  and  Figure  12 
shows  that  the  dominant  flow  is  in  the  plane  perpendicular  to  the  tank  axis. 
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In  conclusion,  a description  of  the  COCHISE  facility  has  been  provided.  Of 
primary  interest  was  the  determination  of  the  pressure  on  the  tank  centerline  as 
a function  of  the  inlet  temperature  and  mass  flow  rate.  The  inlet  temperature  and 
pressures  are  measured  at  a location  in  the  inlet  tube  approximately  25  cm  up- 
stream of  the  tank.  Figure  4 was  constructed  to  determine  the  conditions  under 
which  these  measurements  were  meaningful  representations  of  the  properties  of 
the  gas  where  it  flows  into  the  tank.  Gas  temperatures  could  increase  as  much  as 
20  percent  if  the  constraints  of  Figure  4 are  not  obeyed.  An  analysis  of  the  flow 
on  the  tank  centerline  illustrates  that  the  centerline  pressure  is  within  10  percent 
of  being  uniform  along  the  axis  of  the  tank.  Shock  waves  do  exist  within  the  tank 
and  come  within  1.  7 cm  of  the  tank  centerline.  This  distance  increases  to  above 
2.  5 cm  when  the  gas  inlet  temperature  is  raised  above  200°K.  The  predictions 
for  the  centerline  pressure  are  given  in  Figure  10.  For  any  inlet  temperature  T^, 
there  exists  a maximum  mass  flow  per  inlet  tube  for  which  the  flow  in  the  inlet 
tube  is  laminar.  Higher  mass  flows  would  introduce  turbulent  pipe  flow,  and  its 
effect  has  not  been  assessed.  There  exists  a minimum  mass  flow,  below  which 
diffusion  from  the  chamber  walls  may  affect  the  tank  centerline.  This  effect  need 
not  necessarily  be  detrimental  to  the  facility.  The  constraint  merely  indicates  the 
onset  of  the  phenomenon  and  a detailed  analysis  would  be  required  to  examine  it 
further.  Such  a model  is  discussed  in  Appendix  A. 

The  present  results  are  applicable  only  for  the  case  of  a "cold  wall."  Raising 
the  wall  temperature  would  allow  the  condensed  gas  to  vaporize.  The  maximum 
wall  temperature  for  which  the  present  analysis  is  valid  is  also  shown  in  Figure  10. 
Within  these  constraints,  the  present  analysis  predicts  the  chamber  pressure  as  a 
fimction  of  mass  flow  and  gas  temperature.  A possible  source  of  error  in  the  pre- 
diction will  lie  in  the  "static  gas"  approximation  made  in  Section  2.  5.  From  the 
mass  flow  rates,  we  determine  that  the  Mach  number  at  the  wall  is  approximately 
1/3.  This  has  a direct  effect  on  pressure  of  only  10  percent.  The  cooling  of  the 
gas  in  the  vicinity  of  the  wall  should  have  a direct  effect  on  the  capture  coefficient. 
However,  this  effect  should  be  incorporated  in  Dawson's  data®  of  C as  a function 
of  Tj.  Hence,  using  the  same  model  as  Dawson  should  yield  the  most  accurate 
results.  The  above  prediction  should  yield  10  percent  accuracy  due  to  the  static 
incompressible  assumption.  Since  the  incompressible  assumption  was  made  in 
both  regions  3 and  4,  20  percent  errors  may  result. 
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3.  DIAGNOSTIC  AND  DATA  ANALYSIS 


One  of  the  reactions  of  primary  interest  to  the  COCHISE  program  is 
k' 

N('^D)  + Og  NO(v)  + O . (44) 


Reaction  (44)  is  exothermic  by  =3.6  eV,  and  it  is  desired  to  determine  the  fraction 
of  the  reaction  exothermicity  that  is  channeled  to  NO  vibrational  energy  as  well  as 
the  rate  constants  for  creating  NO  in  variovis  vibrational  levels.  The  species 
N(  D)  may  be  created  by  passing  N2  or  a mix  of  N2/Ar  through  a microwave  dis- 
charge. It  is  propKJsed  to  measure  the  N(^D)  concentration  via  absorption  of 

g 

nitrogen  lamp  radiation,  as  described  by  Lin  and  Kaufman.  There  are  three 

2 2 0 

strong  nitrogen -emission  lines  arising  from  N(3  S P)  -►  N(  D ) transitions,  occur- 
ing  at  1492.  62,  1492.  81,  and  1494.  67  A,  that  can  be  used  in  the  absorption  meas- 
urement. Since  optical  filtering  requirements  are  greatly  simplified  if  all  three  of 

2 

these  lines  can  be  monitored,  it  was  desirable  to  specify  the  relative  N(  D)  detec- 
tion efficiency  when  monitoring  all  three  lines  as  compared  to  that  for  monitoring 
just  the  two  closest  lines. 

Such  a determination  involves  the  evaluation  of  the  transmission  of  a Doppler- 
broadened  lamp  line  through  a Doppler-broadened  absorption  medium  at  a lower 
temperature.  A Doppler  lineshape  may  be  described  by  the  formula 


k 

V 


2 


k e 
o 


(45) 


where  k^  is  the  absorption  coefficient  at  frequency  v,  is  the  absorption  coeffi-  ' 
cient  at  line  center,  and  w,  which  contains  the  lineshape  behavior  with  frequency 
variation,  is  defined  by 


(v  - V ) I 

w = — r — — Jtn2  ; 
°D  ’ 


(46) 


bp,  the  Doppler  line  halfwidth,  is  directly  proportional  to  n/T.  The  higher  tem- 
perature emission  lineshape  and  the  resonant  absorption  lineshape  for  a gas  are 
then  related  by 


8. 


Lin,  C.  L. , and  Kaufman, 


F.  (1971)  J.  Chem.  Phys.  55:3760. 
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(47) 


a e 


Because  a > 1,  the  absorption  lineshape  falls  off  more  rapidly  with  increasing 
(v  - Vq)  than  the  emission  lineshape.  The  lamp  line  emission  may  be  represented 


l^=  J C exp  - (Wg)  dw  = J*  C dw 

-00  -« 


where  C is  proportional  to  the  strength  of  the  emitting  line.  The  transmitted  rad- 
iation is  given  by 


00 

/2  2 

C exp  (-wp  exp  (-k^Nf  exp  (-w^))  dw 


where  N is  the  concentration  of  absorbers  (molecules /cm  ) and  t is  the  absorption 
path  length.  Thus  the  fraction  of  emitted  energy  lost  to  resonant  self  absorption 
is  defined  by 


« 3 

r ^ C.  exp  - (w./q-)^  exp  I exp  (-w?))  dw. 


-io  i=l 


. exp  -(w./o)  dw. 


where  w is  for  the  absorption  lineshape,  and  the  summation  is  over  the  three  pre- 
viously  defined  lines  that  are  taken  to  be  non-overlapping.  Equation  (50)  may  be 
represented  as  a series  expansion  for  small  values  of  k^Ni : 


fcV  f y,rf- 


9.  Mitchell,  A.C.G. , and  Zemansky,  M.W.  (1934)  Resonance  Radiation  and 
Excited  Atoms,  Cambridge,  London,  p 323. 


Equation  (51)  has  been  evaluated  as  a function  of  absorber  line  density  for  the 
three  lines  in  question  using  the  line  strengths,  degeneracies,  etc,,  provided  in 
Ref,  8,  The  predicted  absorption  histories  for  the  cases  of  lines  at  1492.  62  and 
1492,81  A only,  and  for  all  three  lines  are  shown  in  Figures  15  and  14,  respec- 
tively. These  predictions  are  nominally  for  a lamp  temperature  of  600°K  and 
absorbing  gas  temperatures  from  40  to  600°K,  but  the  predictions  are  approxi- 
mately valid  for  other  lamp  and  absorber  temperatures  exhibiting  the  same  tem- 
perature ratio  a.  From  a comparison  of  Figures  13  and  14,  it  can  be  seen  that  the 

O 

inclusion  of  the  transition  at  1494.  67  A causes  only  a slight  decrease  in  diagnostic 

sensitivity.  The  relevant  path  length  of  the  COCHISE  facility  is  100  cm;  thus, 

2 9 3 

N(  D)  concentrations  of  10  /cm  could  be  monitored  if  1 percent  changes  in  lamp 
intensity  can  be  detected.  Lastly,  the  variation  of  absorptance  with  absorber  tem- 
perature can  be  seen  to  be  relatively  weak,  the  absorptance  changing  by  less  than 

30  percent  for  a factor  of  15  change  in  temperature. 

2 

Although  the  N(  D)  diagnostic  has  not  yet  been  implemented  on  the  COCHISE 
facility,  a number  of  preliminary  studies  have  been  performed.  The  first  of  these 
involved  a measurement  of  the  infrared  chemiluminescence  arising  from  reaction 
(44).  In  this  experiment  a mixture  of  10  percent  N2/9O  percent  Ar  was  passed 
through  a microwave  discharge  and  then  mixed  in  the  test  chamber  with  oxygen. 
Chamber  conditions  were  estimated  tobeP  = 4X10  ^ torr,  T = 100°K.  A portion 
of  the  spectra  observed  in  this  mix  is  shown  in  Figure  15.  The  dominant  feature 
if  the  data  is  radiation  from  the  NO  fundamental  vibration/rotation  band  whose 
V = 1 -*■  0 band  origin  is  at  5.  33  pm.  There  is  also  some  radiation  from  an  unknown 
source  observed  at  lower  wavelengths,  5.03—5.  23  pm,  as  can  be  seen  in  Figure  15. 

The  first  step  in  the  analysis  of  this  data  is  to  deduce  the  relative  vibrational 
level  populations  of  NO  corresponding  to  this  spectra.  Because  the  vibrational 
distributions  are  non-Boltzmann,  this  step  requires  an  iterative  comparison  be- 
tween the  data  and  an  artificially  generated  spectrum.  This  spectrum  is  synthe- 
sized from  normalized  "basis  functions,"  l^iv)  (where  v is  in  wavenumbers)  which 
are  the  spectral  bands  arising  from  the  individual  vibrational  transitions,  v^  v - 1. 

The  NO  electronic  ground -state  has  one  quantum  of  orbital  angular  momentum, 
t (and  thus  is  a state),  and  also  exhibits  spin  angular  momentum  from  an  unpaired 
electron.  It  has  been  determined^®  that  in  NO  the  coupling  of  the  spin  and  orbital 
angular  momentum  is  very  strong  (Hund's  coupling,  case  a)*^^  and  remains  strong 
for  all  rotational  levels  populated  in  the  present  experiment.  Thus,  two  total 
angular  momentum  (n)  states  are  allowed,  n = 1/2  and  3/2.  Coupling  between  n 


10,  Billingsley,  F.  P.  (1976)  J.  Molecular  Spectroscopy  61:53. 

11.  Herzberg,  G,  (1951)  Molecular  Spectra  and  Molec\Har  Structure  I,  Spectra  of 
Diatomic  Molecules,  2nded.,  D.  Van  Nostrand  Co. , New  York. 
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and  rotation  is  weak,  and,  consequently,  the  Honl-London  linestrength  formulas 
are  good  approximations  for  the  true  rotational  linestrengths.  J,  the  rotational 
quantum  number,  must  be  greater  than  or  equal  to  the  total  angular  momentum,  fl , 
and  thus  J = 0.  5,  1.  5,  2.  5,  ...  for  the  state,  and  J = 1.  5,  2.  5,  3.5,...  for 

are: 


the  ^2^2  state-  The  formulas  for  the  linestrength,  Sj, 


Sj  = (J*  + n)/J' 

= (2J'  + l)n^/J'(J'  + 1)  (52) 

sJ  = (J'  + 1 + n)(J'  + 1 - n)/(J'  + 1)  . 

It  can  be  seen  from  Eq.  (52)  that  the  Q branch  lines  will  dimmish  rapidly  in  inten- 
sity with  increasing  J,  and  that  the  P branch  lines  will  be  more  intense  that  the  R 
branch  lines. 

The  intensity  of  a single  vibration/rotation  transition  is  specified  by;^^ 


J±1 


. 4 „P,  Q,  R . 

he  i/_  S-,  A , N . 
T J v-»v-l  V,  J 


(2  J + 1) 


v-*v-l 


(53) 


where  is  the  Einstein  coefficient,  is  the  band  origin  wavenumber, 

h is  Planck's  constant,  c the  speed  of  light,  and  is  the  frequency  of  the  tran- 
sition, defined  as  the  difference  in  energy  between  the  upper  and  lower  states  of 
the  transition.  These  energies  (in  cm"^)  are  calculated  with  the  formula: 

J " “e^'"  + 1/2)  - ‘Jg  V'" 

using  spectroscopic  values  (in  cm  ^)  of  u)  = 1903.  6,  u x = 13.97,  B = 1.  7042, 

^ « 6 C C “ 

and  a = .0178.  Higher-order  terms  in  v and  J are  negligible.  is  the  rota- 
tional constant  for  vibrational  level  v,  and  is  defined  as 

Bv  = Be  - Og(v  + 1/2)  . (55) 

J in  Eq.  (53)  is  the  population  of  the  radiating  vibration/rotation  state  and  may 
be  rewritten  as 


Nv  J =-^(2J  + 1)  exp  - (J(J  + l)Byhc/kT) 


(56) 
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where  is  the  rotational  partition  function,  k is  the  Boltzmann  constant,  and  T 
is  the  temperature.  is  the  total  number  of  molecules  in  vibrational  state  v. 

Given  Eqs.  (53)  to  (56)  and  the  Einstein  coefficients  from  Billingsley's  work^® 
shown  in  Figure  16,  the  radiant  intensity  of  a given  transition  per  vibrationally 
excited  molecule,  that  is,  specified  for  all  transitions. 

From  this  tabulation,  the  radiation  signature  of  each  vibrational  "basis  function" 
f ^(v)  is  fotind.  Integration  of  the  intensities  over  a system  slit  function  must  be 
performed  prior  to  fitting.  In  the  COCHISE  experiment,  the  luminescence  is 
monitored  by  a grating  spectrometer  that  has  a triangular  slit  function  having  a 
constant  resolution  in  wavelength  space.  Therefore,  in  the  present  analysis  the 
convolution  of  the  intensity  predictions  over  the  slit  function  is  to  be  (performed 
in  wavelength  space. 

Because  the  trcinsition  linewidth  is  much  narrower  than  the  resolution,  the 
slit  function  integral  may  be  replaced  by  a sum,  and  convolution  then  involves  only 
the  calculation  of  the  slit  function's  fractional  transmission  of  each  rotational 
line's  intensity  at  line  center.  The  nominal  system  resolution  was  taken  to  be 
0.013  pm  and  the  rotational  temperature  assumed  to  be  100°K,  in  accordance  with 
experimental  conditions.  The  effect  of  spin  splitting  was  also  included  in  these 
basis  fimctions.  Coupling  of  the  spin  and  orbital  angular  momentum  in  the  molecule 
causes  the  two  total  orbital  angular  momentum  states  to  be  split  by  121  cm'^. 
Similar  transitions  occurring  within  the  two  doublet  states  lie  at  the  same  fre- 
quencies (within  the  resolution  of  the  experiment).  However,  the  ^1^2 
lower  energy  and  is  preferentially  thermally  populated  by  a factor  of  5.  7 at  100°K. 
Both  ^3/2  lines  are  included  in  the  analysis. 

These  radiation  basis  functions  may  be  related  to  sin  observed  intensity  at  any 
wavelength  by  the  relation 


Ip(X)=S  N^5^(X) 


(57) 


where  the  qusmtities  are  proportional  to  the  populations  in  level  v of  the  excited 
NO  molecules.  Given  the  data  I(X1  and  the  functions  5y(X),  defined  over  some  wave- 
length interval,  the  qusmtities  may  be  determined  by  requiring  that  the  square 
of  the  differences  (Ip(X)  - I(X))  be  minimized.  A computer  code  developed  to  per- 
form this  task  takes  the  data  and  basis  functions  as  input  and  provides  the  relative 
vibrational  populations  as  output. 

In  order  to  provide  a rapid  analysis  of  the  data,  a "first  pass"  of  basis  func- 
tions was  developed  under  the  assumptions  of  n = 1 and  integral  values  of  J.  The 
best-fit  computer  spectra  generated  with  these  basis  functions  is  shown  in  Fig- 
ure 17  in  comparison  with  the  data  of  Figure  15.  The  set  of  curves  shown  at  the 
bottom  of  the  spectra  are  the  individual  basis  functions,  summed  to  provide  the 
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predicted  spectra.  These  individual  basis  functions  represent  the  contribution  of 
each  vibrational  level  to  the  total  intensity.  As  can  be  seen,  the  comparison  be- 
tween curve  fit  and  data  is  quite  good,  surprisingly  so  considering  that  only  approx- 
imate basis  functions  were  used.  Note  that  the  contribution  of  the  v = 1 -•  0 tran- 
sition has  been  overestimated,  because  the  data  shows  underlying  spectral  features 
in  this  wavelength  region  that  are  not  due  to  NO. 

More  detailed  basis  functions  were  then  generated  and  the  effect  of  the  spin 
states  was  examined.  A comparison  of  the  data  with  a curve  fit  performed  using 
half-integral  J values,  but  assuming  all  the  molecules  to  have  n = 1/2 
is  shown  in  Figure  18.  The  predicted  Q-branches  are  much  weaker  than  those 
observed  experimentally.  In  addition,  the  synthetic  spectrum  seems  to  have 
greater  resolution  than  the  data.  Consequently,  a curve  fit  in  which  the  resolution 
was  decreased  to  0.02  /urn  was  generated  and  is  shown  in  Figure  19.  As  can  be 
seen,  the  match  with  the  experimental  data  is  much  better  in  this  case,  and  this 
resolution  was  used  in  all  subsequent  calculations.  The  results  of  a similar  cal- 
culation performed  assuming  all  the  population  to  be  in  the  >^3^2  l®vel  are  pre- 
sented in  Figure  20.  The  predicted  Q branches  are  very  prominent,  as  expected 
from  Eq.  (52).  Thus,  even  though  the  n = 1/2  states  are  preferentially  populated 
by  a factor  of  5.  7,  ’''3^2  transitions  significantly  contribute  to  the  spectrum  in  the 
Q branch  regions.  The  fit  of  a synthetic  spectrum,  from  a calculation  including 
thermally  populated  and  ’''3^2  states,  is  cc .npared  to  the  experimental  data  in 

Figure  21.  As  can  be  seen,  the  comparison  between  data  and  curve  fit  is  once 
again  excellent. 

A modified  version  of  the  basis  function-generation  program  was  also  created 
for  use  on  the  PDP8  computer  used  for  COCHISE  data  acquisition.  Due  to  the 
PDPB's  limited  (16  K)  active  memory  space,  the  basis  functions  are  individually 
calculated,  convolved,  normalized,  and  then  integerized  to  minimize  storage  vol- 
ume. The  least -square -fitting  routine  requires  far  too  much  memory  space  for 
use  on  the  PDP8,  so  a new  algorithm  has  been  developed.  The  method  involves 
the  sequential  fitting  of  the  basis  functions  to  the  data  at  values  near  the  peak  in- 
tensities of  the  basis  function  P and  R branches.  At  the  100°K  temperature  under 
which  the  experinient  is  performed,  only  the  first  10  to  15  rotational  levels  of  the 
vibration  are  populated.  Thus,  only  two  vibrational  basis  fvinctions  have  significant 
intensity  at  most  wavelengths.  The  P branch  of  a vibrational  band  occurs  at  almost 
the  same  wavelength  as  the  R branch  of  the  next  higher  vibration,  and  both  bands 
contribute  to  the  emitted  intensity.  No  such  overlap  occurs  for  the  R branch  of  the 
V = 1 transition.  A fit  to  the  emitted  intensity  at  wavelengths  near  the  peak  of  that 
branch  permits  unique  determination  of  the  population  of  that  level.  The  average 
of  2()  values  calculated  near  the  branch  peak  is  taken  as  the  estimate.  This  popu- 
lation is  then  used  in  determining  the  population  of  the  v = 2 level  by  fitting  the  two 
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basis  functions  to  the  data  at  wavelengths  where  the  overlap  of  the  P and  R 
branches  of  the  two  vibrations  is  maximum.  In  this  manner  the  populations  of  the 
vibrational  levels  are  sequentially  determined.  A feed-back  mechanism  to  permit 
modification  of  the  previous  level's  population,  if  the  values  in  the  two  adjacent 
levels  differ  by  too  great  a factor,  is  included  to  damp  oscillation  of  the  solutions. 
After  all  populations  are  estimated,  a complete  spectrum  is  calculated  using  the 
basis  functions.  A fit  to  the  data  with  8 percent  standard  deviation  has  been 
achieved  using  this  technique.  Hand-fitting  of  the  spectrum,  an  option,  with  slight 
modification  of  the  sequentially  calculated  populations,  has  resulted  in  a 5 percent 
standard  deviation  fit. 

In  conclusion,  two  techniques  have  been  developed  for  determining  NO  vibra- 
tional populations  from  fundamental  band  radiation  signatures.  The  first  of  these 
is  a completely  self  consistent  computer  analysis  that  minimizes  the  square  of  the 
differences  between  data  and  prediction.  This  technique  should  be  used  in  any 
final  data  analysis.  The  second  technique  involves  a computer-interactive  hand- 
fitting technique  that  can  be  used  on  site  during  data  acquisition  to  provide  a rela- 
tively accurate  estimate  of  the  relevant  vibrational  populations. 

The  desired  end  result  of  the  exper- 
iment is  to  relate  the  deduced  vibra- 


tional populations  to  the  rate  constants 
for  reaction  (44).  The  relative  vibra- 
tional populations,  as  determined  from 
the  data  fit  shown  in  Figure  21,  are 
shown  in  Figure  22.  As  mentioned 
earlier,  the  population  of  the  v = 1 state 
has  been  overestimated,  since  the  data 
exhibits  an  additional  radiation  source 
in  the  wavelength  region  of  the  v = l-*0 
transition.  Note,  with  the  exception  of 
level  1,  that  the  vibrational  population 
varies  smoothly,  decreasing  monotoni- 
cally  with  increasing  vibrational  level. 
(The  deduced  populations  for  levels 
V = 12-13  are  only  approximate,  given 
the  large  background  radiation  in  that 
wavelength  region. ) It  is  interesting  to 
note  that  the  relative  vibrational  popula- 


NO VIBRATIONAL  LEVEL 


Figure  22.  Relative  NO  Vibrational 
Populations  as  Deduced  from  the 
Data  of  Figure  15 


tions  deduced  from  the  fits  shown  in  Figures  17  to  20  do  not  differ  significantly 


from  those  shown  in  Figure  22. 
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A kinetic  analysis  is  required  in  order  to  relate  these  vibrational  populations 
to  the  rate  constants  for  reaction  (44).  The  first  step  in  this  analysis  is  to  esti- 
mate the  relative  concentrations  of  the  various  species  in  the  test  chamber.  In  a 

microwave  discharge  the  nitrogen  atom  concentration  will  typically  be  1 percent 

2 2 

of  that  of  N2,  and  the  N(  D),  N(  P)  concentrations  would  be  1 percent  of  the  nitro- 
gen atom  concentration.  Thus,  given  the  total  pressure  of  4 mtorr,  one  can  esti- 
mate the  species  concentrations;  these  are  listed  in  Table  2.  The  large  helium 
concentration  shown  in  Table  2 is  the  result  of  a helium  leak  in  the  system, 
present  at  the  time  of  the  measurement.  The  residence  time  within  the  test 
chamber  is  estimated  to  be  1 msec  (this  estimate  could  be  a lower  bound  because 
of  the  helium  leak);  thus,  chemical  reactions  that  occur  on  a time  scale  of  <10  msec 
will  be  of  interest. 


I 

1. 

► 


1 


Table  2.  COCHISE  Species  Concentrations, 
P = 4 mtorr,  T = lOOOK 


Species 

3 

Concentration,  Part  /cm 

O2 

4 X 10^^ 

Ar 

4 X 10^^ 

^2 

1 2 

3 X 10^^ 

N(^S) 

3 X 10^° 

N(^D,  ^P) 

3 X 10® 

He 

14 

4 X 10  ^ 

The  three  reactions  that  can  produce  NO  in  the  system  are: 


N(^S)  + Og  - NO  + 0 

kgg  = 1.  1 X 10'^^  T exp  (-3150/T)  cm® /sec 

(58) 

N(^D)  + O2  - NO  + 0 

k..  = 5.2  ± .4  X 10  cm®/sec 

44  ' 

(44) 

N(^P)  + O2  NO  + 0 

kgg  = 2.  6 ± . 2 X 10  cm®/sec  , 

(59) 

12 

where  the  rate  constant  for  reaction  (58)  was  taken  from  Garvin  and  Hampson, 

and  those  for  reactions  (44)  and  (59)  are  room  temperatures  measurements  reported 


12,  Garvin,  D.,  and  Hampson,  R.F.,  eds.  (1974)  Chemical  Kinetics  Data  Survey 
VII  Tables  of  Rate  and  Photochemical  Data  for  Modeling  of  the  Strato- 
sphere (Revised),  NBS1R74-430,  January. 
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13 

by  Husain  et  al.  Reactions  (44)  and  (59)  are  also  the  only  important  loss  mech- 
2 2 

anisms  for  N(  D,  P).  Reaction  (58)  may  be  neglected  at  the  temperatures  of 

interest,  given  its  high  activation  energy.  On  the  other  hand,  both  reactions  (44) 

and  (59)  may  have  contributed  to  the  observed  NO  spectra. 

The  peak  concentrations  of  NO  and  O that  could  build  up  within  the  gas  resid- 

7 3 

ence  time  are  estimated  to  be  6 X 10  /cm  : thus,  any  reactions  between  NO(v)  and 
NO  or  O will  be  unimportant.  The  remaining  reactions  that  can  deactivate  vibra- 
tionally  excited  NO  are: 


NO(v)  NO(v  - 1)  + hv 

(60)  ! 

NO(v)  + N Ng  + O k = 2.2  X 10“^^  cm^/sec 

(61) 

NO(v)  + N2  -»•  NO(v  - 1)  + Ng 

(62) 

NO(v)  + O2  ->  NO(v  - 1)  + O2 

(63) 

NO(v)  + Ar  -»•  NO(v  - 1)  + Ar 

(64) 

NO(v)  + He  ^ NO(v  - 1)  + He  . 

(65) 

The  Einstein  coefficients  for  NO(v)  are  shown  in  Figure  16,  and  it  can  be  seen  that 
radiative  decay  will  be  unimportant  since  > 15  msec  for  all  vibrational 

levels.  Furthermore,  the  characteristic  time  for  reaction  (61)  is  approximately 
1 sec,  thus  this  reaction  may  be  neglected.  The  vibration  to  translation  (or  vibra- 
tion to  vibration)  reactions  (62)  to  (65)  are  more  difficult  to  assess  inasmuch  as  the 
rate  constants  for  these  processes  increase  with  vibrational  level  and  are  for  the 
most  part  not  well  defined,  particularly  for  T = 100°K.  Shown  in  Table  3 is  the 
available  rate  constant  information  for  these  processes  as  well  as  the  rate  constant 
that  would  be  required  to  provide  a characteristic  reaction  time  of  10  msec  for  the 
species  concentrations  listed  in  Table  2.  As  can  be  seen,  it  appears  that  reactions 
(62)  to  (65)  will  not  cause  vibrational  deactivation  of  NO  during  the  gas  residence 
t time. 

In  conclusion  then,  vibrationally  excited  NO  will  be  created  by  reactions  (44) 
and  (59)  and  the  only  important  loss  mechanisms  will  be  fluid  dynamical.  There- 
fore, the  relative  NO  vibrational  populations  shown  in  Figure  22  are  directly  pro- 
I portional  to  the  NO(v)  creation  rates.  The  relative  contributions  of  N(  D)  and 

i 

13.  Husain,  D. , Miera,  S.  K. , and  Young,  A.N,  (1974)  J.  Chem.  Soc.  Far. 

Trans.  II  70:1721. 
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Table  3.  Rate  Constant  Data  for  Vibrational  Deactivation  of  NO 


1 


Process 

Necessary  Rate 
Constant 
k(T  = 10  msec) 
cm^/sec 

Reported  Rate 
Constant 
k cm^/sec, 

T = 300OK 

Ref. 

NO(v)  + He 

2.  5 X 10'^^ 

V = 1 1.  2 X 10'^® 

14 

V = 15  3.4  X 10'^^ 

Estimated 

NO(v)  + Ar 

2.  5 X 10"^^ 

<He 

15 

NO(v)  + Og 

2.  5 X 10'^^ 

V = 1 2.4  X 10'^^ 

16 

V = 7 < 1.4  X 10'^® 

17 

NO(v)  + N2 

3.3  X 10"^^ 

V = 1 1.  5 X 10'^® 

14 

V = 1 1.  7 X 10'^® 

16 

2 2 2 
N(  P)  to  these  creation  rates  cannot  be  ascertained  until  such  time  as  the  N(  D,  P) 

absorption  diagnostic  is  implemented. 

Some  effort  has  also  been  directed  towards  determining  the  source  of  the 
underlying  radiation  observed  in  Figure  15  in  the  wavelength  region  of  5.03  to 
5.  23  pm.  This  radiation  exhibits  a band-like  structure  and,  although  not  shown, 
this  structure  extends  down  to  =4.  5 pm  where  an  apparent  bandhead  is  observed 
(the  data  are  quite  noisy  in  this  region  and  this  interpretation  is  open  to  question). 
The  NgOfpj)  bandhead  is  at  4.  5 pm,  thus  it  is  possible  that  the  observed  radiation 
could  be  due  to  highly  vibrationally  excited  NgO,  that  is,  NgO  (OOn).  This  sugges- 
tion remains  to  be  confirmed  by  additional  measurements;  however,  N2O  is  one  of 
the  few  molecules  that  can  be  formed  in  mixtures  of  active  nitrogen  and  oxygen. 

The  N2O  formation  mechanism  remains  to  be  identified.  The  accepted  neutral 
species  reaction  mechanism  for  formation  of  N2O  in  active  N2/O2  mixtures  involves 
a sequence  of  binary  and  tertiary  reactions  and  cannot  be  operative  at  the  densities 
and  time  scales  of  the  COCHISE  experiment.  One  possible  N2O  formation'mech- 
anism  appropriate  to  the  conditions  at  the  COCHISE  facility  involves  the  metastable 


14.  Stephenson,  J.C.  (1973)  J.  Chem.  Phys.  59:1523. 

15.  Wray,  K.  L.  (1962)  J.  Chem.  Phys.  36:2597. 

16.  Murphy,  R.E.,  Lee,  E.T.P,,  and  Hart,  A.M,  (1975)  J.  Chem,  Phys.  63: 

2919.  

17.  Whitson,  M.  E.  Jr.,  Darnton,  L.A.,  and  McNeal,  R.J.  (1976)  Chem.  Ph5rs. 

Lett.  41:552.  
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N2(A  S)  state.  The  NglA)  state  is  formed  in  significant  quantities  in  active  nitro- 
gen and  has  been  found  to  be  rapidly  quenched  by  oxygen  molecules  with  a room 

“’12  3 18 

temperature  rate  constant  of  6 X 10  cm  /sec  . Furthermore,  it  has  been 
18  19 

found  ’ that  oxygen  atoms  are  formed  in  this  quenching  reaction.  One  possible 
branch  for  this  reaction  would  be 

N2(A)  + Og N2O  + O AH=-2.5eV  , (66)  j 

1 

which  is  sufficiently  exothermic  to  allow  the  product  N„0  to  be  created  in  high  5 

vibrational  states.  Reaction  (66),  or  its  equivalent,  will  provide  the  dominant 

-3 

chemical  loss  mechanism  for  N2(A)  and  would  have  a time  constant  of  4 X 10  sec  1 

in  the  test  chamber.  ^ 

The  dominant  formation  mechanism  for  N2(A)  state  is  through  the  nitrogen 
atom  recombination  reaction 

• I 

N + N + M->N2  + M . (67) 

20 

It  is  estimated  that  approximately  one  half  the  nitrogen  molecules  created  by 

reaction  (67)  are  formed  in  the  A state.  Since  this  is  a three-body  reaction,  N2(A)  ' 

state  will  not  be  produced  in  significant  quantities  in  the  test  chamber.  It  will, 

however,  be  formed  in  the  microwave  discharge,  and  its  concentration  in  the  inlet 

pipe  will  be  controlled  by  the  formation  mechanism  (67)  and  quenching  by  the 

reaction 


N2(A)  + N(‘^S)  ->■  N2  + N(^P)  , (68) 

“11  3 8 

which  has  a room  temperature  rate  constant  of  5 X 10  cm  /sec.  Therefore, 
in  steady  state,  the  N2(A)  concentration  may  be  crudely  estimated  to  be  defined  by 

N2(A)  = 0.  5 kg,^  (N)  (M)/kgg  (69) 

where  (N)  and  (M)  are  the  nitrogen  atom  and  total  particle  concentrations,  respec- 

_3 

tively.  The  estimated  inlet  tube  pressure  is  1 to  2 torr;  thus,  for  k__  = 4.4  X 10 

6,  12 

cm  /sec 


18.  Meyer,  J.A.,  Klosterboer,  D. , and  Setser,  D.  W.  (1971)  J.  Chem.  Phys. 

^:2084. 

19.  Meyer,  J.A.,  Setser,  D.  W. , and  Stedman,  D.  H.  (1970)  J.  Phys.  Chem. 

74:2238.  

20.  Thrush,  B.A.,  and  Wild,  A.H.  (1972)  J.  Chem.  Soc.  Far.  Trans.  I 68:2023. 
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N2(A)  = 3 X 10“®  W'^S)  . (70) 

Since  this  ratio  would  be  maintained  in  the  expansion  from  the  inlet  tube  to  the 
center  of  the  test  chamber,  the  estimated  N,(A)  concentration  in  the  chamber 

CO  ^ 

would  be  =10  part/cm  . Therefore,  over  the  gas  residence  time  of  1 msec, 

2,4  X 10^  part/cm^  of  N2O  could  be  produced  by  reaction  (66).  This  is  to  be  com- 
pared to  the  concentration  of  NO  produced  by  reactions  (44)  and  (59),  which  has 
been  shown  to  be  6 X 10  /cm  . While  the  estimated  N2O  concentration  is  more 
than  3 orders  of  magnitude  below  that  for  NO,  it  is  to  be  noted  that  NgOfvg)  rad- 
iates 25  times  more  strongly  than  NO,  whereas  the  unidentified  band  radiation  is 
less  than  one  tenth  that  of  the  NO  radiation. 

The  crude  order  of  magnitude  arguments  presented  here  indicate  that  it  is 
possible  that  N2O,  produced  by  reaction  (66),  could  be  the  cause  of  the  underlying 
band  radiation  seen  in  Figure  15.  This  conjecture  must  be  tested  by  more  detailed 
measurements  with  more  definitive  data  analysis.  Note  that  reaction  (66)  is  of 
considerable  interest  in  its  own  right  inasmuch  as  it  would  provide  an  upper 
atmospheric,  or  auroral,  source  for  N2O  production. 

Two  other  preliminary  measurements  have  been  performed  on  the  COCHISE 
facility  and  these  will  be  discussed  briefly.  In  the  first  of  these  measurements, 
discharge  excited  N2/Ar  was  mixed  in  the  chamber  with  a gas  mixture  composed 
of  10  percent  CO/90  percent  02-  Since  the  discharge  excited  nitrogen  will  be 
highly  vibrationally  excited,  it  was  anticipated  that  the  reaction 


N2(v)  + CO  - N2(v  - 1)  + CO(l) 


would  occur  and  the  resulting  CO  fluorescence  could  be  used  to  provide  an  internal 
wavelength  calibration  for  the  system.  The  observed  fluorescence  in  this  case  is 
shown  in  Figure  23.  The  peak  to  the  left  of  the  figure  corresponds  to  the  v - 1 -►  0 
transition  of  CO,  whose  band  center  is  at  4,  666  ;um.  Thus  it  was  found  that  the 
nominal  wavelength  calibration,  shown  in  Figure  23,  was  off  by  0.  033  pm.  This 
correction  was  used  in  NO  data  analysis. 

More  interestingly,  it  can  be  noted  that  CO  radiation  from  the  transitions 
v=2-».ltov=5-4  was  also  observed.  These  states  could  have  been  excited  by 
multi-quantum  vibrational  transfer,  that  is, 


N2(v)  + CO  -»  N2(v  - n)  + CO(n) 
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or  alternately  by  a sequence  of  single  quantum  exchange  reactions  such  as 


N2(v)  + CO(v')  -»  N2(v  - 1)  + CO(v'  + 1) 


(73) 


and 


CO(v)  + CO(v')  ^ CO(v  - 1)  + CO(v'  + 1)  . (74) 

Given  the  low  density  and  short  residence  time  within  the  test  chamber,  it  is  antic- 
ipated that  reaction  sequences  involving  multiple  collisions  between  species  pro- 
duced only  in  the  test  chamber  are  unlikely,  and  that  the  observed  fluorescence  is 
more  probably  due  to  multi-quantum  vibrational  transfer.  Thus  it  is  expected  that 
additional  measurements  of  this  type,  taken  under  more  carefully  controlled  con- 
ditions, will  provide  estimates  of  the  relative  probabilities  of  single  and  multiple 
quantum  transfer  collisions  in  CO. 

The  last  measurement  to  be  discussed  was  again  made  for  calibration  purposes. 
In  this  case  discharge -excited  N2/Ar  was  mixed  with  Ar  in  the  test  chamber  in 
order  to  determine  background  radiation  levels  in  the  absence  of  02*  The  rather 
surprising  spectrum  observed  in  this  case  is  shown  in  Figure  24.  After  a review 
of  possible  radiation  sources,  it  was  determined  that  this  spectrum  is  due  to  fluo- 
rescence from  the  fundamental  vibration-rotation  band  of  CN  and  exhibits  radiation 
from  vibration  levels  as  high  as  v = 12.  Thus,  a carbon-containing  species  must 
have  been  present  as  an  impurity  in  the  gas  mix.  Furthermore,  this  impurity 
must  react  with  active  nitrogen  in  such  a manner  as  to  directly  produce  highly 
vibrationally  excited  CN.  Although  this  impurity  cannot  be  uniquely  identified  at 
present,  there  is  a possibility  that  it  was  CO  inasmuch  as  CO  had  been  one  of  the 
test  gases  in  the  run  performed  prior  to  this.  No  CO  vibro -luminescence  was 
seen,  however,  and  it  is  pHDSsible  that  the  CN  was  produced  from  another  species, 
for  example  carbonaceous  impurities  in  the  microwave  discharge.  The  CN  pre- 
cursor can  be  determined  in  future  measurements. 

21  -24 

CN  formation  has  been  previously  observed  in  mixtures  of  active  nitro- 

gen and  CO.  In  these  cases,  however,  the  CN  "red"  (A^ir  -*  X ^)  amd  "violet" 

(B^  -»  X ^S)  systems  were  monitored  rather  than  the  infrared  emission.  The 

production  of  CN  in  mixtures  of  active  nitrogen/CO  is  of  fundamental  interest  to 

23 

areas  such  as  laser  physics,  and  several  conflicting  mechanisms  have  been 

21.  Broida,  H.  P. , and  Heath,  D,  F.  (1957)  J.  Chem.  Phys.  ^:1352. 

22.  Stair,  A.T,  Jr.,  Kennealy,  J.P.,  and  Murphy,  R.  E.  (1967)  J.  Chem.  Phys. 

2^:52. 

23.  Taieb,  G. , and  Le gay,  F.  (1970)  Can.  J.  Phys.  48:1956. 

24.  Young,  R.A.,  and  Morrow,  W.  (1974)  J.  Chem.  Phys.  60:1005. 


the 

(75) 

CN(X  ^ CN(A  + N2  (76) 

<=)  CN(B  ^S)  + Ng 

to  be  operative  in  discharge  excited  mixtures  of  Ng/CO.  The  quantity  in 

reaction  (76)  represents  highly  vibrationally  excited  N„. 

24  ^ 

On  the  other  hand  Young  and  Morrow  studied  CN  formation  in  mixtures  of 
active  nitrogen  and  CO  and  concluded  that  the  relevant  CN  formation  mechanism 
was 

N2(B^’r  or  ^A^)  + CO -»  CN(X^S)  + NO  (77) 

followed  by 

N2(^S)  + CN(X^S)  ^ N2  + CN(B^S)  . (78) 

25 

Washida  et  al.  rejected  this  mechanism,  arguing  that  the  rate  constant  for  reac- 

“193  26 

tion  (77)  must  be  «3  X 10  cm'^/sec.  Furthermore,  Provencher  and  McKenney 

have  independently  concluded  that  reaction  (76),  rather  than  (78),  provides  the  dom- 

24-26 

inant  excitation  mechanism  for  CN.  The  only  point  on  which  all  authors  agree 
is  that  N2(A)  state  is  not  the  source  of  CN  excitation. 

It  is  suggested  that  measurements  on  the  COCHISE  facility  would  be  particu- 
larly valuable  in  evaluating  the  various  proposed  CN  formation  and  excitation 
mechanisms.  The  basic  advantage  of  the  COCHISE  facility  is  that  the  densities  are 
so  low  that  only  vibrationally  excited  nitrogen  and  metastable  states  such  as  N2(A) 
will  be  present  in  the  gas  mix.  Therefore,  the  kinetic  interpretation  of  any  ob- 
served fluorescence  is  greatly  simplified. 

For  example,  the  observations  of  Figure  23  have  been  analyzed  under  the 
assumption  that  the  system  impurity  was  indeed  CO  (this  remains  to  be  shown 
experimentally).  The  reactive  species  present  in  the  test  chamber  are  limited  to 

25.  Washida,  N.,  Kley,  D.,  Becker,  K.  H. , and  Groth,  W.  (1975)  J.  Chem.  Phys. 

M:4230.  

26.  Provencher,  G.  M. , and  McKenney,  D.J.  (1972)  Can.  J.  Chem.  50;2527. 
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4 2 2 

N(  S).  N(  D),  N(  P),  N2(A),  and  N2(v).  The  possible  reactions  between  nitrogen 


atoms  emd  CO  are: 

N(^S)  + CO  -►  CN(X  ^S)  + 0 

AH  = 3.3  eV 

(79) 

N(^D)  + CO  CN(X  ^L)  + 0 

AH  = 0.  9 eV 

(80) 

N(^P)  + CO  -►  CN(X  ^S)  + 0 . 

AH  = -0.  3 eV  . 

(81) 

The  first  two  of  these  may  be  discarded  inasmuch  as  they  are  endothermic.  Al- 

13 

though  reaction  (81)  is  exothermic,  it  has  been  shown  that  there  is  no  adiabatic 

2 

reaction  path  to  the  specified  reactants,  (The  quenching  of  N(  P)  by  CO  has,  how- 

13 

ever,  been  found  to  be  rapid,  and  it  is  possible  that  the  reaction  could  proceed 
by  a non-adiabatic  transition.  In  any  event,  the  exothermicity  of  reaction  (81)  is 
insufficient  to  provide  the  CN  vibrational  excitation  observed. ) The  reaction 

N2(v)  + CO -►  CN  + NO  (82) 

requires  5.  5 eV  in  vibrational  energy  to  be  exothermic  and  may  be  neglected. 

18 

Lastly,  the  reaction  between  N2(A)  state  and  CO  has  been  found  to  produce 
CO(a^Jf),  that  is, 

N2(A^S,  V = 0,  1)  + CO  - CO(a^)r)  + N2  k = 2.  2 X 10'“  cm^/sec  . (83) 

Reaction  (83)  is  almost  thermoneutral  and,  thus,  at  room  temperature  the  product 
3 

CO(a  Jr)  state  may  be  rapidly  quenched  by  N2,  again  reforming  N2(A).  This  path 
would  become  much  less  probable  at  the  COCHISE  temperature  of  100°K. 

It  is  concluded  then  that  the  only  species  produced  in  reaction  between  the  dis- 

3 

charge  species  and  CO  is  CO(a  n).  A possible  subsequent  path  for  formation  of 
CN  would  then  be 

N(^S)  + CO(a^ff) -*  CN(X^ff)  + O AH  = -2.  7 eV  . (84) 

It  is  interesting  to  note  that  this  reaction  is  sufficiently  exothermic  to  produce  CN 
vibrationally  excited  to  v a 12,  as  can  be  observed  in  the  data  of  Figure  24.  This 
result  may  well  be  a coincidence,  but,  given  the  present  analysis  as  well  as  the 
controversy  in  the  literature,  it  would  be  of  value  to  attempt  to  duplicate  the  obser- 
vations of  Figure  24  under  controlled  conditions. 


i 


;}! 


4.  SUMMARY  AND  CONCLUSIONS 

In  summary,  a model  has  been  developed  to  describe  the  fluid  d5rnamical  fields 
within  the  COCHISE  facility.  This  model  has  been  exercised  to  predict  expected 
test  chamber  conditions,  such  as  pressure,  temperature,  and  gas  residence  time, 
as  a function  of  system  parameters  such  as  gas  flow  rate  and  wall  temperature.  ^ 

Furthermore,  limiting  system  operating  conditions,  such  as  the  gas  flow  rate  at 
which  the  inlet  pipe  flow  will  become  turbulent,  have  been  designated.  Lastly, 
the  computational  procedure  required  for  a more  detailed  flow  analysis  of  the  sys- 
tem has  been  outlined  in  Appendix  A. 

Analyses  and  interpretation  of  preliminary  infrared  fluorescence  measure- 
ments taken  on  the  system  have  also  been  presented.  The  most  prominent  of 
these  was  the  analysis  of  NO  fundamental  vibration  band  radiation  resulting  from 
the  chemi -excitation  reaction 


I 

i 


N(^D)  + Og  ->  NO(v)  + O . (44) 

Two  computer  codes  were  developed  that  allow  determination  of  the  NO  vibrational 
populations  corresponding  to  the  observed  spectrum.  The  first  of  these  provides  a f 

best  least -squares -fit  between  the  data  and  a computer  generated  spectrum.  The  f 

second,  developed  for  use  on  the  COCHISE  facility  PDP8  computer,  allows  use  of  { 

a less  accurate,  computer  interactive,  spectrum-fitting  procedure,  (Although 

developed  for  NO,  these  codes  may  be  used  for  fitting  the  measured  fundamental  | 

vibrational  band  spectrum  of  any  diatomic  molecule.)  The  NO  vibrational  popula-  j 

tions  corresponding  to  the  observed  spectrum  were  determined  from  the  first  of  i 

these  codes,  and  it  was  demonstrated  that  these  relative  vibrational  populations  I 

were  directly  proportional  to  the  creation  rates  for  NO(v).  ^ 

Several  other  pieces  of  preliminary  data  were  analyzed  in  terms  of  CO  vibro- 
excitation  and  in  terms  of  unique  formation  mechanisms  for  N2O  and  CN.  This 
data  analysis  was  complicated  both  by  a lack  of  diagnostic  information  and  by  insuf- 
ficient low  temperature  kinetic  data  and,  for  the  most  part,  must  be  considered  as 

pure  conjecture.  Nonetheless,  a number  of  intriguing  spectral  features  have  been  j 

observed  in  this  preliminary  data.  It  is  suggested  that  the  examination  of  such  > 

data,  taken  under  controlled  conditions,  could  provide  a wealth  of  useful  kinetic  ) 

information.  I 

i 


44 


1.  Albini,  F.A.  (1965)  AIAA  Journal  3:1535.  j 

2.  Hubbard,  E.  W.  (1966)  AIAA  Journal  4:1877.  | 

3.  Boynton,  F.  P.  (1967)  ALAA  Journal  5:1703.  ! 

4.  Ashkenas,  H. , and  Sherman,  F.  S.  (1966)  The  structure  and  utilization  of 

supersonic  free  jets  in  low  density  wind  tunnels,  in  Rarefied  Gas  Dynamics,  1 

Fourth  Symposium,  Vol.  II.  Academic  Press,  New  York,  pp.  84-105. 

5.  Bier,  K. , and  Schmidt,  B.  (1961)  Zeitschrift  Agnew  Physik  1^:495.  j 

6.  Dawson,  J.P.  (1966)  J.  Spacecraft  3:218.  j 

7.  Schlichting,  H.  (1955)  Boundary  Layer  Theory,  McGraw  Hill.  i 

8.  Lin,  C.L. , and  Kaufman,  F.  (1971)  J.  Chem.  Phys.  55:3760. 

I 9.  Mitchell,  A.C.G.,  and  Zemansky,  M.W.  (1934)  Resonance  Radiation  and  , 

j Excited  Atoms,  Cambridge,  London,  p 323. 

' 10.  Billingsley,  F.  P.  (1976)  J.  Molecular  Spectroscopy  61:53. 

11.  Herzberg,  G.  (1951)  Molecular  Spectra  and  Molecular  Structure  I.  Spectra  of 

Diatomic  Molecules,  2nd  ed. , D.  Van  Nostrand  Co.,  New  York.  | 

12.  Garvin,  D. , and  Hampson,  R.F.,  eds.  (1974)  Chemical  Kinetics  Data  Survey  | 

Vri.  Tables  of  Rate  and  Photochemical  Data  for  Slodeling  of  the  Strato-  | 

sphere  (Revised),  NBS1R74-430,  January.  ! 

13.  Husain,  D. , Miera,  S.  K. , and  Young,  A.N.  (1974)  J.  Chem.  Soc.  Far.  | 

Trans.  11  70:1721.  | 

14.  Stephenson,  J.C.  (1973)  J.  Chem.  Phys.  ^:1523.  j 

15.  Wray,  K.  L.  (1962)J.  Chem.  Phys.  36:2597. 

r 16.  Murphy,  R.E.,  Lee,  E.T.P.,  and  Hart,  A.M.  (1975)  J.  Chem.  Phys.  63:  , 

' 2919.  ' 

I 17.  Whitson,  M.  E.  Jr.,  Darnton,  L.A.,  and  McNeal,  R.J.  (1976)  Chem.  Phys. 

i:  Lett.  41_;552.  , 

! 


45 


74:2238.  

20.  Thrush,  B.A.,  and  Wild,  A.H.  (1972)  J.  Chem.  Soc.  Far.  Trans.  I 68:2023. 

21.  Broida,  H.  P. . and  Heath,  D.  F,  (1957)  J.  Chem.  Phys.  26:1352. 

22.  Stair,  A.T.  Jr,,  Kennealy,  J.P.,  and  Murphy,  R.E.  (1967)  J.  Chem.  Phys, 

J^:52.  

23.  Taieb,  G, . and  Le gay,  F.  (1970)  Can.  J.  Phys.  48:1956. 

24.  Young,  R.A.,  and  Morrow,  W.  (1974)  J.  Chem.  Phys.  60:1005. 

25.  Washida,  N. , Kley,  D. , Becker,  K.  H. , and  Groth,  W.  (1975)  J.  Chem.  Phys. 

M:4230.  

26.  Provencher,  G.  M.,  and  McKenney,  D.J.  (1972)  Can.  J.  Chem.  50:2527. 


46 


I 


I 

I 

I 


Appendix  A 

Radiol  Flow  Model  for  COCHISE 

The  analysis  of  Section  2 was  conducted  in  order  to  make  predictions  of  the 
tank  centerline  pressure  as  a function  of  wall  temperature,  mass  flow/inlet  tube, 
and  inlet  temperature.  In  the  process  of  conducting  the  analysis,  it  was  assumed 
that  region  four  was  a stationary  gas.  This  approximation  was  made  in  order  that 
the  model  be  compatible  with  the  manner  in  which  the  capture  coefficients  were 
measured.  A more  detailed  model  for  region  four  was  originally  initiated.  The 
model  included  a spatially  variable  mean  gas  velocity,  gas  temperature,  and  gas 
pressure.  When  it  was  realized  that  the  capture  coefficients  were  not  known  accu- 
rately enough  to  justify  such  a model,  the  effort  was  dropped.  The  model  is  docu- 
mented in  this  appendix  for  two  reasons:  first,  inaccuracies  in  the  present  model 
may  force  us  to  look  more  closely  at  the  pressure  and  temperature  variations  in 
the  tank  and,  secondly,  the  possibility  of  gas  diffusing  from  the  chamber  wall  to 
the  tank  centerline  could  be  assessed  more  accurately  using  the  radial  flow  model 
that  we  will  now  describe. 

In  Section  2.  6,  we  determined  that  most  of  the  flow  is  in  the  plane  perpendicu- 
lar to  the  tank  axis.  Figure  12  of  the  main  text  illustrates  that  the  flow  approaches 
the  wall  as  though  it  were  coming  from  a line  source  on  the  tank  centerline.  Hence, 
we  define  a radial  flow  model  with  axial  symmetry.  The  mass  originates  at  a line 
source  on  the  axis  of  the  tank  and  flows  radially  toward  the  sink  at  the  cryogenically 
cooled  walls. 
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Let  us  define  the  reference  temperature,  pressure,  density,  and  velocity  as 


1 


Tref=Ti 

(A-1) 

Pref  = Mu^^f/2xRY 

(A -2) 

Pref  " Pref/^'^i 

(A -3) 

and 

"ref= 

(A -4) 

where  6^  is  the  gas  constant  and  R is  the  radius  of  the  tank. 

The  dimensionless  variables  are  defines  as 

e = T/T. 

(A -5) 

= «/^ref 

(A -6) 

“ = P/Pref 

(A -7) 

^ = p/p  ref 

(A -8) 

and 

a = r/R  . 

(A-9) 

The  Reynolds  and  Knudsen  numbers  become 

Re  = M/p  TT  Y = 0(10) 

(A-10) 

and 

Kn*v/R^=0(^)  , 

(A-11) 

respectively.  The  approximate  values  of  Re  and  Kn  suggests  that  an  appropriate 
description  of  the  tank  flow  can  be  made  with  the  use  of  the  Navier  Stokes  equations 
and  "slip  flow"  boundary  conditions.  ^ 

1.  Schaaf,  S.  A.  (1958)  Flow  of  rarefied  gases,  in  Fundamentals  of  Gas  Dynamics, 
ed.  H.  W.  Emmons,  Princeton  Univ.  Press,  1958,  pp.  709-719. 
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Navier  Stokes  Equations 

The  equation  of  continuity  for  radial  flow  with  axial  symmetry  is  expressed  as 


V = 1 


The  momentum  equations  is  expressed 


^ 1 ^ + i 9v 

^ Ba  2 da  Re  1 g^2  a da 


(A-12) 


V 

T 

O' 


(A-13) 


where  we  have  assumed  that  the  viscosity  is  constant.  The  energy  equation 
becomes 


^ da  \ y / ^ PrRe 


9 d . ±96, 

g^2  a da 


4(7  - 1) 
7 Re 


(1^)  ■"T 


7J  ’ 


(A-14) 


where  the  second  and  third  terms  on  the  right  are  the  conduction  and  dissipation, 
respectively,  and  Pr  is  the  Prandtl  number.  The  equation  of  state  is  expressed 
as 


ij  = <t>6 


(A-15) 


The  density  ^ and  the  pressure  u may  be  eliminated  via  the  continuity  and  state 
equations,  respectively.  Second-order  equations  for  the  velocity  v and  tempera- 
ture 6 result: 


9^v  _v  19v,r.ri9v  0 /9v,\,l  9e'| 

^ ■ a 9^  + [a  ^ r + V ^ 9^J 


(A-16) 


and 


9^0 


da 


^ _ 4Pr(7  - 1) 
a da  7 


. .2  2 

9v\  v 


+ Pr  Re  - 

Q OQ 


-PrR.  (ULl)  i |»  . » 

'7/0'  oa  Q 'V  ao-  o/ 


(A-17) 
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Wall  Boundary  Conditions 


I 

I 

t 

} 


Let  us  assume  that  the  gas  at  the  wall  has  a Maxwellian  velocity  distribution 
characteristic  of  the  gas  temperature  T^,  pressure  p^,  mean  velocity  u^,  and 
density  p^; 


JVl 

i 


(2^6iTJ 


3/2 


exp 


2 2 2 


2(R  T„ 


(A-18) 


where  is  the  molecular  velocity  normal  to  the  wall  and  1 2 ?3  tangent  to 

the  walls.  Let  us  also  assume  that  all  molecules  leaving  the  wall  are  thermally 
accommodated  to  the  wall  temperature  and  obey  a wall  Maxwellian  distribution 
function  characteristic  of  that  temperature; 


f = N exp 
w ^ 


2<5^"w 


(A-19) 


where  N is  to  be  determined. 

The  flux  of  mass  per  unit  area  toward  the  wall  (M^)  is  determined  in  a manner 
described  by  Schaaf  : 

00  00  00 

f^’  d?j  d?2  dCg  (A -20) 

-«  *“00  O 


or 


M 


D F, 
w 1 


(2 


(A-21) 


where 

Fj  = exp  (-/3^)  +v/>  p 3 erf /3  . (A-22) 

$ is  the  speed  ratio, 

3 = u /{2(R.T  , (A-23) 

' w ' w ' 
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^ 2 

erf  p = J e * dx  , 
'/'X  J 


(A -24) 


In  a similar  manner,  we  obtain  the  flux  of  translational  energy  toward  the  wall: 


00  00  00 

/ If  ^l<?'/2)f^d|^d€2dC3 


(A -25) 


-00  - 00  O 


M A / A'  <■» 

tr  w 2 


(A -26) 


where 


F„  = (1  + |3^/2)  exp  (-p^)  + (3^  + 5/2))3(l  + erf  p).J'v/2  . 


In  addition  to  translational  energy,  each  molecule  carries  j internal  modes: 


j = (5  - 3y)/(7  - 1)  . 


(A-27) 


(A -28) 


Each  mode  possesses  internal  energy  5^T^/2.  Hence,  the  flux  of  internal  energy 


^int  = ^ ^ ■ 


(A -29) 


The  mass  flux  away  from  the  wall  is  denoted  by  M • 


00  00  o 

iff 


(A-30) 


*00  - 00  - 00 


M‘  = 2:rN((RT^)^  . 


(A-31) 
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The  translational  energy  flux  from  the  wall  is 


00  00  o 


tr 


- f f f 


(A-32) 


"00  - 00  - 00 


(A-33) 


and  the  flux  of  internal  modes  is 


E 


- =j(^) 

int  V 2 / 


int 


M 


(A-34) 


From  the  positive  and  negative  mass  fluxes,  we  define  the  sticking  coefficient 
S in  much  the  same  way  that  the  capture  coefficient  is  defined,  except  that  the  mass 
flux  terms  now  contain  the  mean  flow  u : 


S = 


M^  - M‘ 


M 


(A -35) 


From  the  mass  and  energy  fluxes  at  the  wall,  we  express  the  conservation  of  mass 
and  energy  through  the  tank.  Hence, 


(M  -M)2JrRY  = 2M 


(A-3G) 


and 


Kr^<nt-Et‘r-^rnt)2’^«Y  = 2J^CpT.  , 


(A-37) 


where  T.  is  approximately  the  total  temperature  of  the  inlet  gas. 

Equations  (A-35)  through  (A-37)  are  three  independent  equations  in  N,  T 


and  p^.  Their  solutions  yield 


^w/’^i  = ^3  <1  - + 27SF3/(y  + 1) 


(A -38) 


and 


Pw/Pref=2vA’^^/T;7T-/SF^  , 


(A-39) 


i 

ll 


I 


where 


F3  = (4  + j)/(4F2/F^  + j)  . (A-40) 

The  equations  of  state  and  continuity  are  used  to  relate  the  speed  ratio  to  the 
sticking  coefficient, 

S = 2N/'7r  |3/Fj  . (A-41) 

Equations  (A -38)  through  (A -41)  express  the  boundary  conditions  to  be 
imposed  on  the  gas  at  the  wall.  The  gas  temperature  T^,  pressure  p^.  and 
velocity  u^  (or  speed  ratio  ^)  are  expressed  as  functions  of  the  wall  temperature 
T , inlet  gas  temperature  T.,  sticking  coefficient  S,  and  the  mass  flow  per  inlet 
tube  M. 

Before  proceeding  with  the  analysis,  let  us  examine  the  wall  boundary  condi- 
tions in  more  detail.  Figure  A-1  illustrates  the  gas  temperature  as  a function  of 
sticking  coefficient  for  specific  values  of  wall  temperature.  For  values  of  S near 
zero,  the  gas  assumes  the  wall  temperature.  This  result  v/as  built  into  the  model 
when  it  was  assumed  that  the  molecules  that  don't  stick  to  the  wall  are  thermally 
accommodated.  As  the  sticking  coefficient  increases,  molecules  are  removed 
from  the  tank  and  not  replaced  with  cold  molecules  from  the  wall.  Hence,  the  gas 
temperature  increases  with  sticking  coefficient.  This  trend  continues  until  the 
sticking  coefficient  becomes  sufficiently  close  to  unity  that  the  wall  offers  no  back 
pressure.  The  gas  expands  as  though  it  were  expanding  to  a vacuum,  and  the  gas 
temperature  decreases  in  a supersonic  expansion.  The  speed  ratio  P approaches 
infinity  as  the  sticking  coefficient  approaches  unity.  Figure  A -2  illustrates  the 
gas  pressure  as  a function  of  gas  temperature  T^,  wall  temperature  and 
sticking  coefficient  S.  As  S ♦ 0,  the  pressure  increases  beyond  limit  because  we 
assumed  a constant  influx  of  mass  M.  For  increasing  values  of  S,  the  pressure 
decreases  and  temperature  increases  until  we  approach  S > 0.  9,  where  the  gas 
expands  supersonically  toward  the  wall  with  decreasing  temperature  and  pressure. 

Source  Boundary  Conditions 

In  order  to  integrate  Eqs.  (A-Ui)  and  (A-17)  from  the  tank  centerline  to  the 

wall,  we  must  specify  four  boundary  conditions  on  v and  0.  We  have  determined 

V and  0 and  must  now  determine  two  conditions  on,  or  near,  the  tank  centerline, 
w w 

If  we  consider  the  flux  of  mass  to  be  transonic  at  some  r-adius  o , then  the  bound- 
ary conditions  become 


I 


0 .1  .2  .3  .4  .5  .6  .7 

S - Sticking  Coefficient 


Figure  A-1.  Gas  Temperature  at  Wall 


T n.  ■ Gas  Temperature  at  Wall 


Figure  A -2.  Gas  Pressure  at  Wall 
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and 


I 

1 

I 


i 


0 = e*  = 2/(7  + 1) 


(A -43) 


at 


(A -44) 


where  is  the  unknown  value  of  Po3/Pref’  quantity  that  we  are  attempting  to 
calculate.  Mathematically,  we  can  consider  adjusting  until  one  wall  boundary 
condition  is  satisfied.  Integration  of  the  governing  equations  requires  the  specifi- 
cation of  the  gradients  in  0 and  v at  o = a . Since  the  Navier  Stokes  equations 
were  introduced  to  describe  the  heat  conduction  at  the  wall,  we  assume  that  the 
transonic  flow  at  the  source  is  inviscid  and  the  gradients  are  related  by 


-2(7  - 1) 
y 


* 9 V 

da 


a 


* 


(A -45) 


Since,  in  the  inviscid  limit,  both  gradients  are  infinite  at  the  sonic  point,  we  must 
choose  some  bounded  value  of  dv/da  to  start  the  integrations  and  adjust  that  grad- 
ient until  the  second  wall  boundary  condition  is  satisfied. 


Computational  Procedure 

1)  Specify  7,  Pr,  M and  p (Re  = M/p’rY). 

2)  Specify  Arbitrary  Starting  Conditions. 

i)  Uq  - Pog/Pref 

ii)  grad*  = dv/da  at  u*. 

3)  Determine  Initial  Conditions. 

i)  V*  from  (A  -42) 

ii)  0*  from  (A  -43) 

iii)  O’*  from  (A -44) 

iv)  90/9o  from  (A-45). 

4)  Integrate  (A -16)  and  (A -17)  from  o = o*  to  o = 1.  Obtain  v^  and  0^  at 
0=1. 

5)  From  the  definition  of  3,  3 - v /\^0  . 

' ' w'  w 

6)  Calculate  erf  3,  F^,  F2  and  Fg  as  defined. 

7)  Determine  .S  from  (A -41), 


I 


I 


I 
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8) 


Determine  r 


compatible  with  (A -3  8): 


T.  (1  - S)F, 


2 7 =1 

(y  + 1)(1  - S) 


9)  Repeat  for  an  array  of  values  of  grad. 

10)  Parametrically  replace  grad*  with  the  wall  temperature  t^/T.. 

11)  Have  generated,  at  fixed  M,  y,  Pr,  andp: 


ref 


S 


Fct 


Fct 


Fct 


12)  Require  data  or  theory  for  S(t^,  T^,  u^). 

13)  Steps  11)  and  12)  ^4  relationships  for  S,  u^,  and  in  terms  of  the 

known  values  of  T.  and  t . Hence  we  now  have  Poo^'^o  t,..,  M). 

1 w 2 3 ^ ” 

Unfortunately,  step  12  is  a bottleneck.  Dawson's  data  gives  S(T.,  and 

I w 

that  data  is  insufficient  for  this  model.  More  detailed  models  for  the  sticking 
coefficients  are  required  before  the  radial  flow  model  is  justified. 
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LIBRARY  BRANCH  G-043 
WASHINGTON,  DC  20545 
ATTN  DOC  CON  FOR  CLASS  TECH  LIB 

DEPARTMENT  OF  TRANSPORTATION 
OFFICE  OF  THE  SECRETARY 
TAO-44,1,  ROOM  10402-R 
400  7TH  STREET  S.W. 

WASHINGTON,  DC  20590 
ATTN  SAMUEL  C CORONITI 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MO  20771 
ATTN  A C AIKEN 

NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
ATTN  A TEMPKIN 


NASA  ' 
600  INDEPENDENCE  AVENUE  S W | 
WASHINGTON,  DC  20546  ! 
ATTN  R FELLOWS  | 

NASA  ! 
600  INDEPENDENCE  AVENUE  S W ■. 
WASHINGTON,  DC  20546  | 
ATTN  A SCHARDT  | 

NASA  I 
600  INDEPENDENCE  AVENUE  S W j 
WASHINGTON,  DC  20546  I 
ATTN  M TEPPER  i 

NASA  i 
LANGLEY  RESEARCH  CENTER  I 
LANGLEY  STATION  i 
HAMPTON,  VA  23365  j 
ATTN  CHARLES  SCHEXNAYDER  MS-168 
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NASA 

AMES  RESCH  CENTER 
MOFFETT  FIELD,  CA  90435 
ATTN  N- 254-4  WALTER  L.  STARR 

NASA 

AMES  RESEARCH  CENTER 
MOFFETT  FIELD.  CA  94035 
ATTN  N-254-4  R WHITTEN 

NASA 

AMES  RESEARCH  CENTER 
MOFFETT  FIELD,  CA  94035 
ATTN  N-254-4  ILIA  G POPPOFF 

NASA 

AMES  RESEARCH  CENTER 
MOFFETT  FIELD,  CA  94036 
ATTN  N-254-3  NEIL  H FARLOW 

NASA 

GEORGE  C MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  AL  35812 
ATTN  C R BALCHER 

NASA 

GEORGE  C MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  AL  35812 
ATTN  H STONE 

NASA 

GEORGE  C MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  AL  35812 
ATTN  W A ORAN 

NASA 

GEORGE  C MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  AL  35812 
ATTN  CODE  ES22  JOHN  WATTS 

NASA 

GEORGE  C MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  AL  35812 
ATTN  W T ROBERTS 


NASA 

GEORGE  C MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE,  AL  35812 
ATTN  R D HUDSON 

NASA 

GEORGE  C MARSHALL 
SPACE  FLIGHT  CENTER 
HUNTSVILLE.  AL  35812 
ATTN  R CHAPPELL 

ALBANY  METALLURGY  RESCH  CENTER 
U S BUREAU  OF  MINES 
P.O.  BOX  70 
ALBANY.  OR  97321 
ATTN  ELEANOR  ARSHIRE 

CENTRAL  INTELLIGENCE  AGENCY 
ATTN  RD/SI  RM  5G48  HQ  BLDG 
WASHINGTON  DC  20505 
ATTN  NED/OS I -2G4R  HQS 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  DC  20234 
ATTN  SEC  OFFICER  FOR 
ATTN  JAMES  DEVOE 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  DC  20234 
ATTN  SEC  OFFICER 
STANLEY  ARRAMOWITZ 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON.  DC  20234 
ATTN  SEC  OFFICER  FOR 
ATTN  J COOPER 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  DC  20234 
AHN  SEC  OFFICER  FOR 
AHN  GEORGE  A SINNATT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  DC  20234 
AHN  SEC  OFFICER  FOR 
ATTN  K KESSLER 
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DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON.  DC  20234 
ATTN  SEC  OFFICER  FOR 
ATTN  M KRAUSS 


AERODYNE  RESEARCH.  INC 
BEDFORD  RESEARCH  PARK 
CROSBY  DRIVE 

BEDFORD.  MA  01 731  ATTN  M CAMAC 


DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON.  DC  20234 
ATTN  SEC  OFFICER  FOR 
ATTN  LEWIS  H GEVANTMAN 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON.  DC  20234 
ATTN  SEC  OFFICER  FOR 
ATTN  JAMES  DEVOE 

NATIONAL  OCEANIC  & ATMOSPHERIC  ADMIN 

ENVIRONMENTAL  RESEARCH  LABORATORIES 

DEPARTMENT  OF  COMMERCE 

BOULDER.  CO  80302 

ATTN  GEORGE  C REID  AERONOMY  LAB 

NAT  OCEANIC  & ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  FO  COMMERCE 
BOULDER.  CO  80302 
AHN  ELDON  FERGUSON 


AERONOMY  CORPORATION 
217  S NEIL  STREET 
CHAMPAIGN.  IL  61820 
ATTN  A BOWHILL 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES.  CA  90009 
ATTN  N COHEN 

AEROSPACE  CORPARATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  HARRIS  MAYER 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  SIDNEY  W KASH 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  T WIDHOPH 
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NAT  OCEANIC  & ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  FO  COMMERCE 
BOULDER.  CO  80302 
ATTN  FRED  FEHSENFELD 

AERO-CHEM  RESCH  LABS,  INC 
P.O.  BOX  12 
PRINCETON,  NJ  08540 
ATTN  A FONTIJN 

AERO-CHEM  RESCH  LABS.  INC 
P.O.  BOX  12 
PARINCETON.  NJ  08540 
ATTN  H PERGAMENT 

AERODYNE  RESEARCH,  INC 
BEDFORD  RESEARCH  PARK 
CROSBY  DRIVE 

BEDFORD,  MA  01731  ATTN  F BIEN 


AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  R J MCNEAL 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  R GROVE 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  IRVING  M GARFUNKEL 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
AHN  THOMAS  D TAYLOR 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES.  CA  90009 
ATTN  V JOSEPHSON 
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AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  JULIAN  REINMEIMER 


BROWN  ENGINEERING  COMPANY  INC 
COMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 
ATTN  RONALD  PATRICK 


AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
ATTN  R D RAWCLIFFE 

AVCO-EVERETT  RESCH  LAB  INC 
2385  REVERE  BEACH  PARKWAY 
EVERETT,  MA  02149 
ATTN  TECHNICAL  LIBRARY 

AVCO-EVERETT  RESCH  LAB  INC 
2385  REVERE  BEACH  PARKWAY 
EVERETT,  MA  02149 
ATTN  GEORGE  SUTTON 

AVCO-EVERETT  RESCH  LAB  INC 
2385  REVERE  BEACH  PARKWAY 
EVERETT,  MA  02149 
ATTN  C W VON  ROSENBERG  JR 


CALIFORNIA  AT  RIVERSIDE,  UNIV  OF 
RIVERSIDE,  CA  92502 
ATTN  ALAN  C LLOYD 

CALIFORNIA  AT  RIVERSIDE,  UNIV  OF 
RIVERSIDE,  CA  92502 
ATTN  JAMES  N PITTS  JR 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
3175  MIRAMAR  ROAD 
LA  JOLLA,  CA  92037 
ATTN  S C LIN 

CALIFORNIA  UNIVERSITY  OF 
BERKELEY  CAMPUS  ROOM  318 
SPROUL  HALL 
BERKELEY,  CA  94720 

ATTN  SEC  OFFICER  FOR  HAROLD  JOHNSTON 


BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 
COLUMBUS,  OH  43201 
ATTN  DONALD  J HAMMAN 


CALIFORNIA  UNIVERSITY  OF 

BERKELEY  CAMPUS  ROOM  318 

SPROUL  HALL 

BERKELEY,  CA  94720 

ATTN  SEC  OFFICER  FOR  F MOZER 


BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 
COLUMBUS,  OH  43201 
ATTN  DONALD  J HAM 

BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 
COLUMBUS,  OH  43201 
ATTN  STOIAC 


BATTELLE  MEMORIAL  INSTITUTE 
505  KING  AVENUE 
COLUMBUS,  OH  43201 
ATTN  RICHARD  K THATCHER 

BROWN  ENGINEERING  COMPANY  INC 
CUMMINGS  RESCH  PARK 
HUNTSVILLE,  AL  35807 
ATTN  N PASSING 


THE  TRUSTEES  OF  BOSTON  COLLEGE 
CHESTNUT  HILL  CAMPUS 
CHESTNUT  HILL,  MA  02167 
ATTN  CHAIRMAN  DEPT  OF  CHEM 


CALIFORNIA  UNIVERSITY  OF 

BERKELEY  CAMPUS  ROOM  318 

SPROUL  HALL 

BERKELEY,  CA  94720 

ATTN  SEC  OFFICER  FOR  DEPT  OF  CHAM 

W H MILLER 

CALIFORNIA,  STATE  OF 
AIR  RESOURCES  BOARD 
9528  TELSTAR  AVENUE 
EL  MONTE,  CA  91731 
ATTN  LEO  ZAFONTE 

CALSPAN  CORPORATION 
P.O.  BOX  235 
BUFFALO,  NY  14224 
ATTN  C E TREANOR 

CALSPAN  CORPORATION 
P.O.  BOX  235 
BUFFALO,  NY  14221 
ATTN  G C VALLEY 
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CALSPAN  CORPORATION 
P.O.  BOX  235 
BUFFALO.  NY  14221 
ATTN  M G DUNN 

CALSPAN  CORPORATION 
P.O.  BOX  235 
BUFFALO,  NY  14221 
ATTN  W WURSTER 

COLORADO,  UNIVERSITY  OF 
OFFICE  OF  CONTRACTS  AND  GRANTS 
380  ADMINISTRATIVE  ANNEX 
BOULDER,  CO  80302 
ATTN  A PHELPS  JILA 

COLORADO,  UNIVERSITY  OF 
OFFICE  OF  CONTRANCTS  AND  GRANTS 
380  ADMINISTRATIVE  ANNEX 
BOULDER,  CO  80302 
AHN  JEFFREY  B PEARCE  LASP 

COLORADO,  UNIVERSITY  OF 
OFFICE  OF  CONTRANCT  AND  GRANTS 
380  ADMINISTRATIVE  ANNEX 
BOULDER,  CO  80302 
AHN  C BEATY  JILA 

COLORADO,  UNIVERSITY  OF 
OFFICE  OF  CONTRACTS  AND  GRANTS 
380  ADMINISTRATIVE  ANNEX 
BOULDER,  CO  80302 
ATTN  C LINEBERGER  JILA 

COLORADO,  UNIVERSITY  OF 
OFFICE  OF  CONTRACTS  AND  GRANTS 
380  ADMINISTRATIVE  ANNEX 
BOULDER,  CO  80302 
ATTN  CHARLES  A BARTH  LASP 

COLUMBIA  UNIVERSITY,  THE  TRUSTEES 
IN  THE  CITY  OF  NEW  YORK 
LA  MONT  DOHERTY  GEOLOGICAL 
I OBSERVATORY-TORREY  CLIFF 

PALISADES,  NY  19064 
ATTN  B PHELAN 

COLUMBIA  UNIVERSITY,  THE  TRUSTEES 
OF  THE  CITY  OF  NEW  YORK 
116TH  STREET  & BROADWAY 
NEW  YORK,  NY  10027 
; ATTN  RICHARD  N ZARE 


COLUMBIA  UNIV,  THE  TRUSTEES  OF 

CITY  OF  NEW  YORK 

116TH  & BROADWAY 

NEW  YORK,  NY  10027 

ATTN  SEC  OFFICER  H M FOLEY 

CONCORD  SCIENCES 
P.O.  BOX  119 
CONCORD.  MA  01742 
ATTN  EMMETT  A SUTTON 

DENVER,  UNIVERSITY  OF 
COLORADO  SEMINARY 
DENVER  RESEARCH  INSTITUTE 
P.O.  BOX  10127  DENVER,  CO  80210 
ATTN  SEC  OFFICER  FOR  MR  VAN  ZYL 

DENVER,  UNIVERSITY  OF 

COLORADO  SEMINARY 

DENVER  RESEARCH  INSTITUTE 

P.O.  BOX  10127  DENVER,  CO  80210 

ATTN  SEC  OFFICER  FOR  DAVID  MURCRAY 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
ATTN  DASAIC 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
ATTN  WARREN  S KNAPP 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER) 
SANTA  BARBARA,  CA  93102 
ATTN  TIM  STEPHENS 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
ATTN  DON  CHANDLER 

GENERAL  ELECTRIC  COMPANY 
TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
ATTN  B CAMBILL 
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GENERAL  ELEC.  CO. 

SPACE  DIVISION 

VALLEY  FORGE  SPACE  CTR 

GODDARD  BLVD 

KING  OF  PRUSSIA 

P.O.  BOX  8555 

PHILADELPHIA,  PA  I9I01 

ATTN  M H BORTNER,  SPACE  SCIENCE  LAB 

GENERAL  ELEC.  CO. 

SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD.  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
ATTN  J BURNS 

GENERAL  ELEC.  CO. 

SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  LOX  8555 
PHILADELPHIA,  PA  19101 
ATTN  F ALYEA 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
ATTN  D HENDERSON 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
ATTN  J S WAGNER  PHYSICS  DEPT 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
ATTN  B J WATKINS 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
ATTN  T N DAVIS 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
ATTN  R PARTHASARATHY 


GENERAL  ELEC.  CO. 

SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
ATTN  P ZAVITSANDS 

GENERLA  ELEC.  CO. 

SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
ATTN  R H EDSALL 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
ATTN  NEAL  BROWN 

LOWELL,  UNIVERSITY  OF 

CENTER  FOR  ATMOSPHERIC  RESEARCH 

450  AIKEN  STREET 

LOWELL,  MA  01854 

AHN  G T BEST 

LOCKHEED  MISSILES  & SPACE  COMPANY 

3251  HANOVER  STREET 

PALA  ALTO,  CA  94394 

ATTN  JOHN  KUMER  DEPT  52-54 


GENERAL  ELEC.  CO. 

SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
ATTN  T BAURER 

GENERAL  RESEARCH  CORPORATION 
P.O.  BOX  3587 
SANTA  BARBARA,  CA  93105 
ATTN  JOHN  ISE  JR 


LOCKHEED  MISSILES  & SPACE  COMPANY 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  KIMER  DEPT  52-54 

LOCKHEED  MISSILES  & SPACE  COMPANY 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  JOHN  B CLADIS  DEPT  52-12 
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LOCKHEED  MISSILES  & SPACE  CO 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  TOM  JAMES  DEPT  52-54 

LOCKHEED  MISSILES  & SPACE  CO 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  J B REAGAN  D/52-12 

LOCKHEED  MISSILES  & SPACE  CO 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  MARTIN  WALT  DEPT  52-10 

LOCKHEED  MISSILES  & SPACE  CO 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  RICHARD  G JOHNSON  DEPT  52-12 

LOCKHEED  MISSILES  & SPACE  CO 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  ROBERT  D SEARS  DEPT  52-14 

LOCKHEED  MISSILES  & SPACE  CO 
3251  HANOVER  STREET 
PALO,  ALTO,  CA  94304 
ATTN  J R WINKLER 

INSTITUTE  FOR  DEFENSE  ANALYSE 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
ATTN  ERNEST  BAUER 

INSTITUTE  FOR  DEFENSE  ANALYSE 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
AHN  HANS  WOLFHARD 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
ATTN  D ARCHER 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
ATTN  D FISCHER 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
ATTN  M SCHEIBE 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
ATTN  D SAPPENFIELD 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
ATTN  D SOWLE 

PHOTOMETRIC,  INC.  j 

442  MARETT  ROAD  j 

LEXINGTON,  MA  02173 
ATTN  IRVING  L KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.O.  BOX  1069 
BERKELEY,  CA  94701 
ATTN  J B WORKMAN 

PHYSICAL  DYNAMICS  INC. 

P.O.  BOX  1069 

BERKELEY,  CA  94701  i 

ATTN  A THOMPSON 

PHYSICAL  SCIENCES,  INC.  ^ 

30  COMMERCE  WAY  i 

WOBURN,  MA  01801 

ATTN  KURT  WRAY  ^ 

3 

PHYSICAL  SCIENCES,  INC.  ] 

30  COMMERCE  WAY  1 

WOBURN,  MA  01801  ! 

ATTN  R L TAYLOR  j 

PHYSICAL  SCIENCES,  INC. 

30  COMMERCE  WAY 

WOBURN,  MA  01801  i 

ATTN  G CALEDONIA 

PHYSICS  INTERNATIONAL  COMPANY 
2700  MERCED  STREET 
SAN  LEANDRO,  CA  94577 
ATTN  DOC  CON  FOR  TECH  LIB 

PITTSBURGH,  UNIV  OF 
OF  THE  COMWLTH  SYS  OF  HIGHER  EOUC 
CATHEDRAL  OF  LEARNING 
PITTSBURGH,  PA  15213 
ATTN  WADE  L FITE 
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PITTSBURGH,  UNIVERSITY  OF 
OF  THE  COMWLTH  SYS  OF  HIGHER  EDUC 
CATHEDRAL  OF  LEARNING 
PITTSBURGH,  PA  15213 
AHN  MANFRED  A BIONDI 

PITTSBURGH,  UNIVERSITY  OF 
OF  THE  COMWLTH  SYS  OF  HIGHER  EDUC 
CATHEDRAL  OF  LEARNING 
PITTSBURGH,  PA  15213 
ATTN  FREDERICK  KAUFMAN 

PITTSBURGH,  UNIVERSITY  OF 
OF  THE  COMWLTH  SYS  OF  HIGHER  EDUC 
CATHEDRAL  OF  LEARNING 
PITTSBURGH,  PA  15213 
ATTN  EDWARD  GERJUOY 

PRINCETON  UNIV,  THE  TRUSTEES  OF 
FORRCSTAL  CAMPUS  LIBRARY 
BOX  710 

PRINCETON  UNIVERSITY 
PRINCETON,  NO  08540 
ATTN  ARNOLD  0 KELLY 

R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
ATTN  RICHARD  LATTER 

R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
AHN  R G LINDGREN 

R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
ATTN  BRYAN  GABBARD 

R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
ATTN  H A DRY 

R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
ATTN  ROBERT  E LELEVIER 

R & D ASSOCIATES 


R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
ATTN  ALBERT  L LATTER 

R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
ATTN  FORREST  GILMORE 

R & D ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
ATTN  D DEE 

R & D ASSOCIATES 
1815  N.  FT.  MYER  DRIVE 
IITH  FLOOR 
ARLINGTON.  VA  22209 
ATTN  HERBERT  J MITCHELL 

R & D ASSOCIATES 
1815  N.  FT.  MYER  DRIVE 
IITH  FLOOR 
ARLINGTON,  VA  22209 
ATTN  J W ROSENGREN 

RAND  CORPORATION 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
ATTN  CULLEN  CRAIN 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 
ATTN  DANIEL  A HAMLIN 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 
AHN  DAVID  SACHS 

SPACE  DATA  CORPORATION 
1331  SOUTH  26TH  STREET 
PHOENIX,  AZ  85034 
ATTN  EDWARD  F ALLEN 

STANFORD  RSCH  INSTITUTE  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK.  CA  94025 
ATTN  M BARON 


P.O.  BOX  9695 

MARINA  DEL  REY.  CA  90291 

ATTN  R P TURCO 
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STANFORD  RSCH  INSTITUTE  INTL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
ATTN  L LEADABRAND 

STANFORD  RSCH  INSTITUTE  INTL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
ATTN  WALTER  H CHESTNUT 

STANFORD  RSCH  INSTITUTE  INTL 
1611  NORTH  KENT  STREET 
ARLINGTON,  VA  22209 
ATTN  WARREN  W BERNING 

STANFORD  RSCH  INSTITUTE  INTL 
1611  NORTH  KENT  STREET 
ARLINGTON,  VA  22209 
ATTN  CHARLES  HULBERT 

TECHNOLOGY  INTL  CORPORATION 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 
ATTN  W P BOQUIST 

UNITED  TECHNOLOGIES  CORP 
755  MAIN  STREET 
HARTFORD,  CT  06103 
ATTN  H MICHELS 


VISIDYNE,  INC. 

19  THIRD  AVENUE 
NORTHWEST  INDUSTRIAL  PARK 
BURLINGTON,  MA  01803 
ATTN  HENRY  J SMITH 

VISIDYNE,  INC. 

19  THIRD  AVENUE 
NORTHWEST  INDUSTRIAL  PARK 
BURLINGTON,  MA  01803 
ATTN  J W CARPENTER 

VISIDYNE,  IND. 

19  THIRD  AVENUE 
NORTHWEST  INDUSTRIAL  PARK 
BURLINGTON,  MA  01803 
ATTN  WILLIAM  REIDY 

VISIDYNE,  INC. 

19  THIRD  AVENUE 
NORTHWEST  INDUSTRIAL  PARK 
BURLINGTON,  MA  01803 
AHN  T C OEGGES 

VISIDYNE,  INC. 

19  THIRD  AVENUE 
NORTHWEST  INDUSTRIAL  PARK 
BURLINGTON,  MA  01803 
ATTN  CHARLES  HUMPHREY 


UNITED  TECHNOLOGIES  CORP 
755  MAIN  STREET 
HARTFORD,  CT  06103 
ATTN  ROBERT  HBULLIS 

UTAH  STATE  UNIVERSITY 
LOGAN,  UT  84321 
ATTN  DORAN  BAKER 

UTAH  STATE  UNIVERSITY 
LOGAN,  UT  84321 
ATTN  KAY  BAKER 

UTAH  STATE  UNIVERSITY 
LOGAN,  UT  84321 
AHN  C WYATT 


WAYNE  STATE  UNIVERSITY 
1064  MACKENZIE  HALL 
DETROIT,  MI  48202 
ATTN  PIETER  K ROL 
CHAM  ENGR  & MAT  SCI 

WAYNE  STATE  UNIVERSITY 
1064  MACKENZIE  HALL 
DETROIT,  MI  48202 
AHN  R H KUMMLER 

WAYNE  STATE  UNIVERSITY 
DEPT.  OF  PHYSICS 
DETROIT,  MI  48202 
AHN  WALTER  E KAUPPILA 


YALE  UNIVERSITY 

UTAH  STATE  UNIVERSITY  NEW  HAVEN,  CT  06520 

LOGAN,  UT  84321  ATTN  ENGINEERING  DEPARTMENT 

AHN  D BURT 
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